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(57) Abstract 

Techniques are provided herein for reducing vibrations 
in various modes of a dynamic system. One such technique 
comprises incorporating vibration limiting and sensitivity 

constraints into a partial fraction expansion equation model 
of the system so as to reduce vibrations to specific levels. 
Another technique comprises shaping a command determined 
using the partial fraction expansion equation model to produce 
a desired output. The entire command may be shaped 
or only selected portions thereof which produce vibrations. 
Another technique involves commanding in current to produce 
saturation in voltage. By doing this, it is possible to command 
voltage switches. The times at which the switches occur 
can be set to reduce system vibrations. Other techniques are 
also provided. These include vaiying transient portions at the 
beginning, middle and/or end of a move and using Posicast 
inputs, among others. 
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DYNAMIC SYSTEM CONTROL METHOD 

BACKGROUND OF THE INVENTION 

Field Of The Invention 

The present invention is directed to altering inputs to, and 
generating feedforward signals for, a dynamic system so as to reduce unwanted 
vibrations in the system. The invention has particular utility in speeding up 
computer disk drives by reducing unwanted vibrations which, if unchecked, could 
lead to disk read/write errors or excessive noise. 

Description Of The Related Art 

Movement in dynamic systems typically results in unwanted 
vibrations that are both mechanical and acoustic in nature. These vibrations can 
have a detrimental affect on the operation of such systems. One dynamic system 
that is particularly sensitive to unwanted vibrations is a computer disk drive. 

A computer disk drive includes an actuator arm having a head 
mounted at a distal end of the arm for reading from, and writing to, tracks on a 
magnetic disk. This head is moved by the arm from track-to-track on the disk. 
Vibrations in the system result from this movement. That is, the head and 
actuator arm vibrate after reaching a nominal final position. Reading and writing 
cannot take place by the disk drive until these vibrations go below a certain level; 
otherwise read/write enors occur. These vibrations may also cause noise. 

The period during which the disk drive waits for these vibrations 
to reach an acceptable level (i.e., the settling time) increases the disk drive's seek 
time. The drive's seek time comprises the lime it takes for the drivers head to 
come to rest at a position where the head can perform a read/write operation on 
a particular track. The increase in seek time can be especially acute in cases 
where the track-per-inch density of the magnetic disk is high. That is, because 
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tracks in these disks are relatively close together, even snniail vibrations in the 
head can seriously affect the accuracy of the disk drive and/or increase noise in 
the drive. Since even small vibrations cannot be tolerated, the settling time is 
further increased, thereby further increasing the drive's seek time. 
5 Conventional attempts at addressing the foregoing problems in both 

disk drives and dynamic systems in general have faUen short of satisfactory. That 
is, such attempts are too computationally intensive to be practical, have failed to 
provide sufficient reduction in vibrations for use in high accuracy positioning 
equipment such as computer disk drives, produce sub-optimal trajectories, and/or 

10 are overly sensitive to system parameter variations. 

Other related problems also plague conventional disk drives. For 
example, in conventional disk drives, two different controllers are used to 
position the drive's head on a track. A first controller controls the drivers head 
to reach a predetermined position near to a final position, at which point a second 

15 controller takes over. This second controller moves the head into the final 
position and regulates the head on a track. Switching between these two 
controllers increases settling time and, as a result, increases seek time. Also, in 
conventional disk drives, little or no control over vibrations is included in the 
derivation of these controllers. Consequently, conventional disk drives may not 

20 be able to discern important variations during motion transients. 

In view of the above, there exists a need for a way to control 
computer disk drives and, more generally, dynamic systems, which reduces both 
mechanical and acoustic vibrations to an acceptable level without undue 
computational effort and without a substantial reduction in movement speed. 

25 

SUMMARY OF THE INVENTION 
The present invention addresses the foregoing needs by providing 
methods of controlling movement in a dynamic system so as to reduce mechanical 
and/or acoustic unwanted vibrations in the system. These methods are relatively 
30 easy to perform, produce good results, and have few adverse side effects, thereby 
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making Ihem a significant improvement over the prior art. 

In one aspect, the invention determines which parameter will cause 
a dynamic system to saturate. For example, in the case of a computer disk drive, 
current commands control the system. However, voltage is the fundamental 
5 limiting parameter of the system, since voltage supplied by the drive's power 
supply limits the amount of current in the system. In this aspect, the fundamental 
limiting parameter is used to generate a command for the system which reduces 
vibrations without causing the system to saturate. In the specific case of a disk 
drive, the fundamental limiting parameter (i.e., voltage) is also used to generate 
10 current commands that do not exceed the voltage limitations of the system. 

In another aspect of the invention, a disk drive system is modeled 
using panial fraction expansion equations. The digital form of those equations 
is as follows: 

N 

Finalpos = V.A^t 

i=\ ' b - a 
,.i ' a - b 

where Finalpos is a final position of the system, A is a scaling factor, a and b are 
15 inverse time constants, tg^d is the time at which a move is completed, Vj are 
yoltage inputs to the system, Tj are the times at which Vj are input, and At is a 
time interval at which are input. Using the above equations, system inputs 
which produce reduced vibrations can be determined. These equations are 
particularly powerful when used to solve for a system input in terms of its 
20 fundamental limiting parameter. Additional constraints may also be included in 
the above equations in order to provide a more accurate model of the system. 

In still anodier aspect of the invention, techniques are provided for 
reducing system vibrations in its various modes. These techniques may be used 
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alone or in conjunction with the concepts set forth above. One such technique 
comprises incorporating vibration limiting and sensitivity constraints into the 
partial fraction expansion equation model so as to reduce vibrations to specific 
levels. Another technique comprises shaping (e.g., convolving or fihering) a 
5 • command determined using the partial fraction expansion equation model to 
produce a desired output. The entire command may be shaped or only selected 
portions thereof which produce vibrations. Selective shaping is preferred since 
it reduces seek time relative to shaping the entire command. Another technique 
involves commanding in current to produce saturation in voltage. By doing this, 

10 it is possible to command voltage switches. The times at which the switches 
occur can then be set to reduce system vibrations. Other techniques are also 
provided. These include varying transient portions at the begiiming, middle 
and/or end of a move and using Posicast inputs, among others. 

In still another aspect, the invention determines an input to the 

15 system that will result in limited (or reduced) vibrations, and uses that input as 
a feedforward trajectory for the system. Several approaches are provided for 
generating the 'ujpnt, including performing an optimization method on the partial 
fraction expansion equations. In one approach, the problem of generatuig an 
input is separated into a rigid mode problem and a flexible mode problem. 

20 Specifically, terms associated with the rigid mode of the system are determined 
from the partial fraction expansion equations, and terms associated with 
oscillating or flexible modes of the system are determined based on a system 
analysis. The rigid body terms are solved for an input which drives the system 
so as to satisfy its rigid body constramts. This input is then shaped using, e.g., 

25 an Input Shaper™ designed to compensate for the flexible modes of the system. 
Using this approach, the computational difficulty of obtaining a solution is 
reduced relative to the optimization approach, while still providing an adequate 
reduction in .levels of vibrations in relevant modes of the system. Another 
approach to generating the input involves creating a closed-loop model of the 

30 system. An input which results in reduced vibrations may then be generated 
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based on this model. In still another approach, voltage switch limes for a bang- 
bang input that reduces vibrations are determined based on the partial fraction 
expansion equations. Input Shaping*™ may then be applied to that input so as to 
provide a further reduction in vibrations. 
5 In other aspects of the invention, various position- velocity ("PV") 

table formulations are provided for reducing system vibrations. In one approach, 
a trajectory which reduces system vibrations is generated in accordance with one 
or more of the techniques set forth above, and that trajectory stored in a PV 
table. The trajectory can then be used to control movement of the system. In 

10 another approach, a shaper is provided at the output of a PV table, and "its inverse 
is provided along the system's feedback path. This arrangement reduces system 
dynamics via the feedforward direction, and restores loop stability in the feedback 
direction. In still another approach, an N (N>2) dimensional PV table is used 
to control a system, in which velocity is a function of position and one or more 

15 other variables, such as movement distance, resistance, temperature, etc. In yet 
another PV approach, current commands are shaped at transition(s) (e.g., at 
deceleration) so as to yield desired levels of vibrations. During non-transitional 
states, a trajectory in a PV table is followed. 

While the above-summarized methods significantly enhance the 

20 operation of dynamic systems, such as computer disk drives and other data 
storage devices, it is also preferable to allow end-users some control over the 
system's operation. For example, a user may want a disk drive that has a 
relatively low seek time (i.e., a disk drive that is fast), but may not care about 
acoustics. Similarly, the user may be more focused on acoustics than on seek 

25 lime. Thus, there exists a need for a way by which an end-user can control 
operation of a disk drive. In panicular, there exists a need for a simple way to 
control operadon of the disk drive, which permits the user to select operational 
parameters, such as seek time and acoustic levels, in addition to other features. 

The present invention addresses these needs by providing a GUI 

30 and/or mechanical switch(es) through which an end-user may alter operational 
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parameters of a disk drive, including, but not limited to, the drive's seek time, 
noise level, and power consumption. By providing these feamres, the invention 
enables the user to customize operation of a disk drive relatively easily without 
requiring significant hardware and/or software modifications. 

This brief summary has been provided so that the nature of the 
invention may be understood quickly. A more complete understanding of the 
invention can be obtained by reference to the following detailed description of the 
preferred embodiments thereof in connection with the attached drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a perspective view of a computer system through which 
the present invention may be implemented. 

Figure 2 shows the architecture of the computer system. 

Figure 3 shows process steps for generating a graphical user 
interface in accordance with the invention. 

Figure 4 to 9 show graphical user interfaces that can be generated 
by the invention. 

Figure lOA is a cut-away view of a disk drive that is controlled by 
the present invention. 

Figure lOB is a block diagram of circuitry for controlling the disk 
drive in accordance with the invention. 

Figure 1 1 is a feedback diagram of a computer disk drive. 

Figure 12 shows a plot of "good" moves versus "bad" moves for 
different movements made by a computer disk drive. 

Figure 13 is a feedback diagram of a prior art disk drive system. 

Figure 14 is a feedback diagram of a disk drive system in 
accordance with the present invention. 

Figure 1 5 is a block diagram showing the relationship between an 
input and an output of a finite impulse response filter shaping technique. 
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Figure 16 is a block diagram showing the relationship between an 
input and an output of a finite impulse response filter inverse shaping technique. 

Figure 17 is a block diagram showing the relationship between an 
input and output of a tapped continuous delay line filter shaping technique. 

Figure 18 is a block diagram showing the relationship between an 
input and an output of a tapped continuous delay line filter inverse shaping 
technique. 

Figure 19 is a block diagram showing the relationship between an 
input and an output of an infinite impulse response filter shaping technique. 

Figure 20 is a block diagram showing the relationship between an 
input and an output of an infmite impulse response filter inverse shaping 
teclmique. 

Figure 21 shows examples of trajectories in a position-velocity 
table used in embodiments of the present invention. 

Figure 22 depicts parameterization of trajectories in a position- 
velocity table in accordance with the present invention. 

Figures 23 A and 23 B show optimal rigid body voltage commands 
for 1 and 3000 track moves, respectively, for a disk drive system. 

Figure 24 shows a digitized input voltage waveform. 

Figure 25 shows the waveform of Figure 24 conceptualized as a 
series of pulses. 

Figure 26 shows each of the pulses of Figure 25 represented as a 
pair of impulses. 

Figure 27 shows a block diagram of a feedback system. 

Figure 28 shows application of a Laplace differentiation operator 
to inputs to the system shown in Figure 27. 

Figure 29 shows scaled versions of the inputs of Figure 28. 

Figure 30 shows a reduction, using block diagram algebra, of the 
inputs shown in Figure 29. 

Figure 31 shows a reduction of the inputs shown in Figure 30. 
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Figure 32 shows a plot of impulses of a three pulse shaper together 
with a sine wave representation of undamped vibration. 

Figure 33 shows a plot of digital pulse spacing for the three pulse 
shaper together with a sine wave representation of undamped vibration. 
5 Figure 34 shows pulse spacing for a three pulse shaper and a sine 

wave representation of vibration. 

Figure 35 shows a voltage pulse response of a dynamic system. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

10 The present invention comprises computer-executable code (i.e., 

computer-executable process steps) for controlling operation of a disk drive or, 
more generally, any dynamic system. For the specific case of a disk drive, the 
invention includes code which generates a GUI, which accepts inputs into that 
GUI, and which outputs commands to a disk drive causing the disk drive to 

15 conform its operation to die GUI inputs. The invention can be used in any type 
of computer/disk drive system, such as a desktop computer, laptop computer, or 
the like. For the sake of brevity, however, the invention will be described in the 
context of a desk-top computer system only. 

Figure 1 shows a representative embodiment of a computer system 

20 1 on which the invention may be implemented. As shown in Figure 1, personal 
computer ("PC") 2 includes network connection 4 for interfacing to a network, 
such as the Internet, an ATM network, or the like, and fax/modem connection 5 
for interfacing with other remote devices such as a digital camera, digital video 
camera, and the like. PC 2 also includes display screen 6 for displaying 

25 information to a user, keyboard 7 for inputting text and user conmiands, mouse 
9 for positioning a cursor on display screen 6 and for inputting user conmiands, 
and disk drive 10 for reading from and writing to floppy disks installed therein. 
PC 2 may also have one or more peripheral devices (not shown) attached thereto 
for inputting text, graphics, images, or the like, and/or for outputting images. 



-8- 



wo 99/45535 PCT/US99/04910 

Figure 2 shows the internal structure of PC 2. As shown in Figure 
2, PC 2 includes memory 11, which comprises a computer-readable medium such 
as a computer hard disk. Memory II stores data 12, applications 14, and 
operating system 15, among other things. In preferred embodiments of the 
5 invention, operating system 15 is a windowing operating system, such as 
Microsoft® Windows98 or Microsoft® Windows NT; although the invention may 
be used with other operating systems as well. Among the applications stored in 
memory 11 is computer code to implement the present invention, namely disk 
drive engine 16. Alternatively, disk drive engine 16 may be stored in a read-only 

10 memory ("ROM"), or the like, together with other operational firmware for disk 
drive 10. In brief, disk drive engine 16 is comprised of computer-executable 
process steps which, inter alia, generate a GUI for controlling operation of disk 
drive 10 and output commands to disk drive 10 causing disk drive 10 to conform 
its operations to settings in the GUI. A more detailed description of disk drive 

15 engine 16 is provided below 

Also included in PC 2 are display interface 17, keyboard interface 
19, mouse interface 20. disk drive interface 21, computer bus 22, RAM 24, 
processor 25, and printer interface 26. Processor 25 comprises a microprocessor 
or the like for executing applications, such those noted above, out of RAM 24. 

20 Such applications may be stored in memory 11 (as noted above) or, alternatively, 
on a floppy disk in disk drive 10. Processor 25 accesses applications (or other 
data) stored on a floppy disk via disk drive interface 21. 

Application execution and other tasks of PC 2 may be initiated 
using keyboard 7 or mouse 9, commands from which are transmitted to processor 

25 25 via keyboard interface 19 and mouse interface 20, respectively. Similarly, 
application execution may be initialed remotely via, e.g., network interface 4 and 
a network interface card (not shown). Output results from applications running 
on PC 2 may be processed by display interface 17 and then output via network 
interface 4 or, alternatively, displayed to a user on display 6. To this end, 

30 display interface 17 preferably comprises a display processor for forming video 
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images based on video data provided by processor 25 over computer bus 22, and 
for outputting those images to display 6. Output results from other applications 
running on PC 2 may be provided to an attached printer via printer interface 26. 

Turning to disk drive engine 16, this application comprises 
5 computer code to control operation of disk drive 10 either according to one of a 
plurality of well*known conventional methods or according to one or more of the 
proprietary methods described herein. Disk drive engine 16 funher includes code 
for generating a GUI for interactively controlling at least one of a seek time, a 
noise level, and a power consumption of a disk drive. The GUI includes a 
10 controller that is operated by a user so as to alter settings in the GUI for one or 
more of these parameters. Thereafter, disk drive engine 16 generates and outputs 
commands to the disk drive causing it to conform its operation in accordance with 
the altered settings in the GUI. 

15 1.0 User Interface 

Figure 3 is a flow chart for controlling operation of a disk drive 
using disk drive engine 16. To begin, step S301 generates a GUI and displays 
that GUI to the user via display 6. Figure 4 shows such a GUI, labeled 27. 
GUI 27 includes a window 70 for selecting a disk drive to be controlled. That 

20 is, although the computer system shown in Figure 1 includes only one disk drive, 
many systems include more than one disk drive. Window 70 provides a way for 
a user to select which disk drive to control. Figure 9 is a close-up view of 
window 70. A disk drive can be selected simply by pointing and clicking on that 
drive in window 70. The selected drive is then highlighted, as shown, 

25 As shown in Figure 4, GUI 27 also includes controller 29. 

Controller 29 is used to alter a seek time and a noise level of disk drive 10. In 
this particular embodiment, controller 29 comprises a sliding bar 30 which moves 
along a continuum 31 between a relatively high noise level/low seek time 
("Quick") end 32 and a relatively low noise level/high seek time ("Quiet") end 

30 34. Disk drive noise level and seek time vary inversely along the continuum. 
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meaning that, as the noise level of the disk drive progressively increases, the seek 
time of die disk drive progressively decreases. Likewise, as the noise level of the 
disk drive progressively decreases, the seek time of the disk drive progressively 
increases. By moving sliding bar 30, it is possible to vary both the seek time and 
5 noise level parameters simultaneously. 

In the preferred embodiment of the invention, sliding bar 30 moves 
between discrete positions on continuum 31, meaning, e.g., between a first 
position corresponding to high noise level/low seek time ("Quick") end 32, a 
second position corresponding to a mid-point between the two ends on continuum 

10 31, and a third position corresponding to low noise level/high seek time ("Quiet") 
end 34. Alternatively, there may be more than three discrete levels or noise and 
seek time may vary continuously between the two ends of continuum 31. 

The level of noise in a disk drive is due to vibrations at particular 
frequencies, such as 6.6 KHz, 6.0 KHz, 5.0 KHz, 4.2 KHz, and 1.3 KHz, On 

15 the other hand, mechanical vibrations, which reduce seek time, result from 
vibrations at different frequencies, such as 3.3 KHz. By moving the sliding bar 
of controller 29 along the continuum, an end-user is effectively determining which 
of the above frequencies should be reduced and/or how much the reduction should 
be. For example, in a case that Input Shaping™ (described below) is used to 

20 reduce vibrations, by moving sliding bar 30, the user effectively selects one or 
more shaping functions for reducmg vibrations at specific frequencies, and causes 
disk drive engine 16 to convolve those functions with the appropriate disk drive 
inputs in order to implement the desired vibration reduction. 

In GUI 27 the noise level and seek time vary inversely along 

25 continuum 31 such that there is an effective trade-off between seek time and noise 
level, meaning that as seek time increases, noise level decreases and vice versa. 
These variations are not necessarily linear; although they may be. As an 
alternative to sliding bar 30, the GUI may include boxes into which discrete noise 
levels and/or .seek times are entered. For example, as shown in Figure 5, the 

30 GUI may include entry box 37 into which a noise level value is entered and 
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display box 39 which displays a corresponding numerical value of the seek time. 
In this case, the GUI preferably also includes up/down arrows 41 which are used 
for selecting the various noise level values. On the other hand, these values may 
be input directly via mouse 9 and keyboard 7. As another alternative, the seek 
5 time may be settable via entry box 37 and the corresponding noise level displayed 
via display box 39. In either case, the numerical value displayed in the display 
box is calculated by disk drive engine 16 based on the value set in the entry box. 

Returning to Figure 4, GUI 27 may also include a second 
controller, such as check box 42. This controller causes the disk drive to enter 

10 a power-saving mode. The power-saving mode comprises a state during which 
the disk drive reduces its power consumption by a predetermined amount. The 
power-saving mode may be set in conjunction with varying the seek tune and the 
noise level. That is, in the power saving mode, the set seek time and noise level 
are achieved using less power than would otherwise be required. This is 

15 achieved, e.g., by convolving one or more predetermined functions with a 
command input to a voice coil motor in the disk drive. 

On the other hand, a reduction in power consumption may be 
achieved via a controller, such as controller 44 shown in Figure 6. Controller 44 
comprises a sliding bar 45 which moves along a continuum 46 between a 

20 relatively high power consumption/low seek time end 47 and a relatively low 
power consumption/high seek time end 49. As was the case above, sliding bar 
45 may move between discrete values along the continuimi. Thus, in this 
embodiment, disk drive power consumption and seek time vary inversely along 
the continuum, meaning that, as the power consumption of the disk drive 

25 progressively increases, the seek time of the disk drive progressively decreases. 
Likewise, as the power consumption of the disk drive progressively decreases, the 
seek time of the disk drive progressively increases. 

By movmg sliding bar 45 it is possible to vary both seek time and 
power consumption simultaneously. For example, varying power consumption 

30 may be achieved by convolving one or more predetermined functions with a 
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command input to a voice coil motor in the disk drive. Selection of the 
predetermined functions is performed based on the settings of the sliding bar. Of 
course, the functions may be set in other ways as well (as is the case throughout 
this application). For example, the functions may be selected from predetennined 
5 trajectories in a position-velocity ("PV") table, as described herein below. 

Figure 7 shows yet another alternative configuration for the GUI. 
As shown in Figure 7, the GUI could include separate controllers 50 and 52 for 
reducing both power consumption and noise level, respectively, as a function of 
seek time. Thus, (i) as noise level increases, seek time decreases, and vice versa, 

10 and (ii) as power consumption increases, seek time decreases, and vice versa. 
The sliding bars shown in Figure 7 preferably have all the characteristics of that 
shown in Figure 4. As shown hi Figure 8, numerical values for noise level and 
power consiunption may be implemented, as in Figure 5. Here, a numerical 
value is calculated for seek time as well based on the noise level and power 

15 consumption settings, and displayed in the GUI as shown. 

Returning to Figure 4, GUI 27 may also include a preview 
controller, shown at 54, to cause the disk drive to operate using parameter(s) set 
by controller 29. This preview controller may comprise a button or the like, 
which a user can click on, and which causes the disk drive to operate in a 

20 "preview" mode. The preview mode may cause the disk drive to seek between 
various tracks on a disk so that the user can experience the noise level and seek 
time settings in the GUI. If these levels are unacceptable, the user may then 
revise the GUI iiqjuts, and preview tho drive's operation based on such newly-set 
inputs. This process can be repeated as often as desired until the user achieves 

25 the desired settings. Preview mode may also be pre-set to execute automatically 
each tune the GUI settmgs are altered. 

GUI 27 may also include one or more display areas 55, as shown 
in Figure 4. These display areas may provide a numerical indication of the noise 
level of the disk drive, its seek time, and/or its power consumption rate. Display 

30 areas such as these are generally most beneficial in cases where parameters are 
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set using a sliding bar» as is the case, for example, in Figures 4, 6 and 7. 
Likewise, a display area may be provided which informs the user of the amount 
of hours that a computer (e.g., a laptop) can operate with an existing power 
supply. To determine this value, disk drive engine 16 monitors the computer's 
power supply directly, or estimates the amount of power left based on the total 
time the computer has been running with the present supply, as well as other 
relevant parameters. 

Remrning to Figure 3, following step S301, in which a GUI is 
displayed, processing proceeds to step S302. In step S302, the user inputs 
desired settings into the GUI. These settings, of course; will depend upon the 
available options in the GUI, examples of which are described above. For 
instance, the user may input some combination of a new seek tune, power 
consumption rate, noise level, etc. Step S302 may be performed locally or 
remotely via network interface 4. That is, die GUI may be displayed on a remote 
networked PC, and the settmgs therefor applied to the local PC. In any event, 
once the new parameter setting are input in step S302, they may be confirmed via 
"OK" button 60, canceled via "Cancel" button 61, or applied to the disk drive 
via "Apply" button 62. The difference between "OK" button 60 and "Apply" 
button 62 is that "OK** button 60 merely stores the new settings in memory (e.g., 
memory 11), whereas "Apply" button 62 applies the settings to the disk drive, 
overriding any previous hardware or software settings. 

Next, in step S303, commands (e.g., electrical signals) are 
generated for controlling disk drive 10 in accordance with the settings set in the 
GUI. Exacfly how these conunands are generated may vary, depending upon the 
way in which the disk drive is controlled. Assuming, for illustration's sake, that 
the disk drive is being controlled via Input Shaping"^, step S303 comprises 
convolving various inputs to the disk drive with predetermined functions which 
are selected based on settings in the GUI in order to produce disk drive inputs 
that achieve the results specified in the GUI. These functions may be stored, 
e.g., in memory 11 or in another memory on the system. Sections 2 to 11 below 
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describe a variety of methods that may be used to control operation of the disk 
drive, any of which may be used in step S303 based on the GUI sening and other 
factors. 

Once the necessary commands have been generated, processing 
5 proceeds to step S304, wherein these commands are provided to control operation 
of the disk drive. Thereafter, a user may again alter the drive's operation simply 
by calling up a GUI using, e.g., mouse 7, and entering desired operational 
parameters. Following step S304, processing ends. 

At this point, it is noted that the GUI is not limited to the 
10 combinations of features shown in the figures. That is, a GUI may be generated 
in step S301 with includes any one or more of the foregoing features, subject to 
certain constraints inherent in the system. Likewise, additional features not 
described herem, but which are well known to those skilled in the an, may also 
be included in the GUI. As an alternative to the GUI, or to supplement the GUI, 
15 the disk drive may include one or more electro-mechanical manual switches for 
controlling all or part of the drive's operation in the manner set forth above. For 
the purposes of the present mvention, a "jumper" is considered to be within the 
definition of electro-mechanical manual switch. 

20 2.0 System Control 

The following control methods are preferably implemented via 
computer-executable process steps in disk drive engme 16 so as to effect high- 
speed moves in a disk drive, such as drive 10. However, it is noted that these 
methods are not limited to use with disk drives. Rather, the methods described 

25 herein can be used to control any type of dynamic system which moves from one 
state to another state. In fact, any unwanted dynamic behavior that can be 
quantified in a mathematical expression can be reduced by the present invention. 
For example, using the present invention, a disk drive in a Redundant Array of 
Independent Drives ("RAID") can be commanded to move in less time than has 

30 heretofore been possible without substantially exciting neighboring drives and 

-15- 




wo 99/45535 PCT/US99/04910 

while simultaneously reducing acoustic vibrations that result in noise. The 
designer simply need select which vibrations or dynamic behavior is troublesome, 
identify constraints and, using the technology presented herein, choose from 
among several approaches for generating optimal or near optimal moves that 
5 reduce those vibrations or dynamic behavior. As used herein, the term 
"vibrations" can refer to mechanical and/or acoustic vibrations that cause noise. 

Figures lOA and lOB show close-up views of disk drive 10. As 
shown, disk drive 10 includes voice coil motor 72 having a rotor (not shown), 
actuator arm 74, sensor 75, and head 76. These components are controlled via 

10 disk drive engine 16, or portions thereof, running on a controller/processor 73. 
This processor may comprise processor 25 set forth above or, alternatively, a 
separate controller dedicated to the disk drive which receives commands from 
processor 25. Head 76 reads to, and writes from, tracks (not shown) on magnetic 
disk 77 installed in drive 10. Actuator arm 74, on which head 76 is mounted, 

15 controls the motion of head 76 between the tracks on magnetic disk 77. Motor 
72 drives actuator arm 74 in accordance with control signals received from 
processor 73. These control signals are generated based» at least in part, on a 
desu"ed position of head 76 and an actual measured position of head 76 provided 
via signal 79. Movement of head 76 from an mitial position to the desired 

20 position occurs within the device's seek time. The seek time for disk drive 10 
comprises the time it takes for head 76 to move from an initial position to rest at 
a position where head 76 can perform a read/write operation on a particular track 
of magnetic disk 77. 

The seek time for computer disk drive 10 can be reduced by 

25 commanding a component in disk drive 10, such as head 76, to move from one 
location to another while satisfying constraints on mechanical vibrations, acoustic 
vibrations, drive mounting vibrations, or any combination of constraints, physical 
limitations, and/or parameter variations. Several different methods for reducing 
such vibrations are described below in detail. One such technique involves 

30 altering system inputs that are specified as a function of time. Another technique 
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involves generating a position-versus-velocity (hereinafter "PV" or "position- 
velocity") table using dynamic elements, e.g., constraints or the like, so as to 
produce PV table trajectories Which result in reduced vibrations. In this regard, 
the terms "iiq)ut" and "trajectory" may refer to either a reference command for 
a controller to follow or a feedforward signal within the controller. Still another 
technique involves modifying the current industry-standard PV table so that an 
improved PV table is constructed and used. Several methods for generating PV 
tables are also provided, together with ways to use an Input Shaper''"'^ in 
combination with a PV table in order to reduce unwanted vibrations. Also 
presented are several modifications to the foregoing embodiments which have 
differing levels of optimality and ease of implementation. 

Processor 73 determines which of these techniques to apply based 
on a variety of factors, such as user inputs to the GUI and the identity of the disk 
drive. For example, time can be traded off against amount of vibration reduction 
and robustness to parameter variation. In a given application, the amount of 
vibration reduction required to make end-of-move conditions acceptable can be 
a function of the movement distance or other parameter. When this is the case, 
it is advantageous to trade-off shaping time against shaper effectiveness (meaning 
amount of vibration reduction) such that (i) vibrations, both mechanical and 
acoustic, are reduced to an acceptable level (not excessively attenuated), and (ii) 
shaping time is less than it would be if the vibration were excessively attenuated. 
When more vibration reduction is required, a longer shaper can be used; when 
less vibration reduction is required a shorter shaper can be used; and when no 
vibration reduction is required, then no shaper (zero length) can be used. In this 
example, shaper selection was from among a discrete collection of shapers. It is 
also possible to use a single adjustable shaper, e.g. one in which the time interval 
between pulses is shortened; resulting in a desired shorter shaping time and 
reduced vibration reduction. 
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3.0 Fundamental Limiting Parameter 

One of the techniques utilized by the present invention is a method 
for commanding systems that saturate in an optimal manner while still reducing 
unwanted dynamics. Since all systems saturate in some manner, a method for 
setting up a command profiling approach is presented herein which first 
determines a fundamental limiting parameter for the system, second applies either 
a command shaping or whole-trajectory-generating technique to that parameter, 
and third schedules the shaped command to be output to the system. 

The fundamental limiting parameter of a system is the parameter 
that causes the system to saturate. In the case of a computer disk drive, current 
commands control the system. However, the fundamental limiting parameter is 
voltage, since voltage supplied by the drive's power supply limits the amount of 
current in the system. By identifying the fundamental limiting parameter of a 
system in this manner, and processing that parameter in accordance with the 
techniques provided herein, it is possible to provide for greater control over 
system movement than has heretofore been possible. 

In this regard, there are several disadvantages to not using the 
fundamental limiting parameter for command profiling. First, there is no a priori 
way to generate a conmiand signal that does not saturate the system. If the 
system is allowed to saturate a signal at any time after any command profiling 
technique is applied, the vibration (or unwanted dynamics) reduction benefit of 
command profiling will be degraded. The amount of degradation depends on how 
much the signal is saturated. Therefore, sub-optimal trajectories often result. As 
examples, optimal rigid body voltage commands for a disk drive system for 1 and 
3000 track moves are shown in Figures 23A and 23B, respectively (see section 
6.3 below for a discussion of rigid vs. flexible system modes). The optimal 
current command that corresponds to this voltage is highly complex because it 
depends on the state of the system (in this case, velocity, due to the back EMF 
of the drive's motor). Any simple current command would necessarily be sub- 
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optimal if the designer is to guarantee that the system is not to saturate. Section 
3.1 below describes how this is achieved. 

For optimal moves, the system should be conunanded with inputs 
that take the system to samration limits without exceeding these limits. The 
5 second disadvantage of not choosing the fundamental limiting parameter is that 
these saturation limits cannot be determined without including the entire system 
dynamics. Therefore, the solution to the problem is significantly more complex. 
Additionally, once vibrational constraints are mcluded in the problem, the solution 
increases yet again in complexity. Once the proper selection of a fundamental 

10 limiting parameter is made, an input that approaches samration without 
unknowingly trying to exceed the limit can be generated. 

Another approach to handling the presence of nonlinearities, such 
as saturation, along the feedforward path is to place a shaping module after the 
nonlinear element. For example, if a signal can saturate along the feedforward 

15 padi, the signal should be saturated before the shaping operation. In the case of 
a digital-to-analog conversion ("DAC") saturation, the output of the feedforward 
calculations might need to be limited to plus or minus 16,000 counts. If a 
shaping technique is to maintain the vibration reducing benefits, the output should 
be limited to 16,000 counts and then shaped. 

20 Still another approach to handling samration of a parameter, while 

maintaining the vibration reduction benefits of the various techniques presented 
herein, is to pre-saturate a signal. One condition under which this would be 
necessary is if the saturating element is outside of the control of the designer 
(e.g., the samration occurs in a hardware amplifier). In this case, it is 

25 impractical to perform shaping after saturation. Under this condition, the signal 
should be pre-samrated and then shaped. Pre-samration involves determining a 
mathematical expression for the condition under which the system will samrate 
during a move, and saturating the signal before it is shaped. 

This approach of placing the shaping operation after any nonlinear 

30 block is applicable to the technique of putting Input Shapers*^^^ and Inverse 
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Shapers™ in a loop, as presented in section 9 below. All of the same rules apply 
in that any saturation or nonlinear operation must occur outside the path between 
the Shaper™ and Inverse Shaper™. If the loop is cut open and all of the dynamic 
elements between the Shaper"^^ and Inverse Shaper'^'^ are extracted, the dynamic 
5 elements must all be linear (or nearly linear); otherwise the effect of unwanted 
dynamic reduction is degraded. Additionally, the pre- saturation technique can be 
used in situations in which saturatmg elements must necessarily be between the 
Shaper"^^ and Inverse Shaper^"^. 

10 3.1 Using The Fundamental Limiting Parameter While Generating 

Feedforward Signals In A Different Parameter 

In a disk drive, current is usually commanded to a 
transconductance amplifier. The fundamental limiting parameter for a disk drive 
IS is voltage. It is, however, possible to generate current conmiands that maintain 
a limit in voltage. This section explains how to generate such commands. 

The first method for generating a current command, that does not 
exceed the voltage limitations of the system is to solve a set of equations directly 
for the conmiand input which, in this case, is the feedforward signal to the 
20 transconductance amplifier. The following equation is used to generate an input 
that takes the system to zero final velocity: 

E^, = o. (1) 

where the A's comprise the amplitudes of the current commands at time interval 
i, and N comprises the last time interval. The following equation determines the 
velocity, v, at each time interval i 

25 where A comprises the current iiq)uts and C,^ comprises a disk-drive-dependent 
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constant. The final desired position/state of the system, is determined as 
follows 

The voltage lunit constraint, V^^, comprises the maximum voltage tiiat is being 
commanded (usually the drive supply voltage). For a disk drive, this is usually 
on the order of 12V. limits the A value by the following equation: 



where Rj is the pulse response of the amplifier and drive system to the voltage 
command. The voltage pulse response may be determined by outputting a unit 
pulse to the amplifier connected to the disk drive and measuring the voltage 
response of the amplifier. The values of Rj are chosen by means of a multiple 
step process outlined below. 

In this regard. Figure 35 shows a voltage pulse response from a 
computer simulation. The smooth curve is the continuous voltage response. The 
"square step" curve represent the system's digital sampling rate (or output rate, 
see section 7.1 below), DT. The first peak occurs Tpej^^ after the initial pulse is 
issued to the system. The first value of R is determined by measuring this peak 
value and setting the first R equal to this peak value. All subsequent values of 
R are determined by samplmg in DT intervals after this initial peak. The step 
curve in the figure shows the actual values for this particular example. In this 
case, a complex voltage curve is reduced to a relatively small number of sample 
points to be used in one of the several solution techniques described herein. 

Alternatively, the pulse response R can be computed as a 
calibration step in the disk drive or generated at design time. All of the required 
information for generating an optunal current command is included in this pulse 
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response (except for slew rate limits discussed below). The amplifier is 
reasonably linear, therefore, samration of the amplifier is easily predicted from 
the pulse response. Also, the pulse response may be altered as a function of a 
measured or estimated parameter change. 

The effect of the back EMF of the drive's voice coil motor on the 
above calculation may be notable. If the voltage pulse response R is calculated 
analytically (or from a model), die voltage is constant and non-zero at steady state 
because the system reaches a finite velocity from a pulse in acceleration. In 
Figure 35, the pulse response achieves a small positive value (which cannot be 
seen on the plotted resolution). Experimentally, friction and flex circuit spring 
effects eliminate this constant velocity and cause the drive's head to stop moving. 
In this case, a constant term in consistent units equal to KbemfCvscale should be 
added to the voltage pulse response R so that it includes the effect of back EMF 
and allows for a more aggressive command. 

An example computer code written in GAMS that may be used to 
solve equations 1 to 4 above for values of A is provided below in Appendix B. 
By solving the set of equations in Appendix B, a tune optimal, unshaped current 
command can be generated. This set of equations can be solved in many different 
ways. For example, the entire set can be solved as a constramed optimization in 
which the goal is to generate the minimum time solution. This may also be 
performed without optimization by solving for a feasible solution (the first step 
of a Simplex solution, for example) and then reducing the length, N, of the 
solution vector until no answer can be found (a binary search may be used here). 
Another method is to determine a function that predicts the length of the solution 
and then, without any iteration, find a feasible solution. One such prediction 
function comprises the switch time equations set forth in Appendix A, described 
below. Another approach is to fit a curves to such a solution. In this case, an 
exact answer is not needed so long as the curve prediction functions overestimate 
the solution lengths and never underestimate the solution lengths. 
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It is also possible to solve for an optimal, unshaped command in 
real-time which is not in the system's fundamental luniting parameter. For 
example, we can solve for cunent (A) by setting a final value of the fundamental 
limiting parameter, here voltage (V), as follows (of course, the invention is not 
limited to use with current and voltage, any parameters can be used). This 
approach is based on the fact that the maximum current command output value 
is related to previous output values by the following relationship: 



i 



AO) = 



5(1) 



10 where R is the voltage pulse response and is the maxnnum voltage on 
acceleration (or any condition under which Ai^.j> A|) and on deceleration (or 
any condition under which A|+i < Aj). Equation 5 assumes diat the back EMF 
of the drive's voice coil motor is already included in R(j). The back EMF term 
may be added to R(j) as set forth above, if needed. Equation 5 can be further 

15 simplified since the R(j) values become insignificant after only a few samples. 
Consequently, an efficient approximation would only use a few samples of the 
R(j) series. If samples are dropped for an efficient approximation, the back EMF 
of the drive's voice coil motor times the number of dropped samples may be 
added to equation 5 to achieve an even more accurate result if that value is 

20 deemed significant. 

A rigid body current command can be generated from equation 5 
in a step-by-step fashion, where each response depends only on previous values. 
First the system is accelerated. Each current step is computed from equation 5. 
Next, the system is decelerated. This deceleration may be to zero for a constant 

25 velocity coast in the middle of the move or it may be a transition to negative 
deceleration and then to a stop. Deceleration is computed from equation 5 by 
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reversing the sign of V^. Under coast conditions, a second deceleration phase 
will be started (again using equation 5 with a negative value for V^. Lastly, 
when deceleration is complete, the transition back to zero current is governed by 
equation 5 with a positive value, once again, for Wj^^^, 

5 The only unknowns in the previous equations are the times at 

which to transition between regions. There are several ways to solve these 
problems. Shaped™ switch times derived in section 6.9 below may be used as 
reasonable estimates for the transition times; although any curve fit or lookup 
table would also suffice. This lookup is especially useful given that on many disk 

10 drives, one current profile is often currentiy used and scaled for a range of seek 
distances. Additionally, using this approach the final position will not be exact 
and a zero velocity must be guaranteed. The problem of zero velocity means that 
an additional constraint must be used in computing the final values. Equation 5 
defines the limit on the current value; therefore, in deceleration, for example, any 

15 value more negative than that value may be used without exceeding the limit. As 
a result, the final values can be chosen so that the sum of the current values 
equals zero (no residual velocity). The distance traveled will be incorrect (but 
very close). The last step is to scale the waveform linearly by the percentage 
error so that the final value is correct. This operation can even be performed as 

20 part of the next step, which is the shaping operation. The shaper can be scaled 
and then convolved with the rigid body current command. Another option is to 
increase the value of N, thus adding extra points, and directly solving for the 
extra A(i) values such that velocity equals zero and position equals "final 
position". 

25 There are many variations to the above approach which may yield ' 

fewer computations. Also any vibration suppression technology can be 
subsequently applied to the rigid body current command. Furthermore, may 
be changed as the calculation progresses to optimize scaling. For example, with 
the knowledge that the current command is to be shaped (and thereby reduced in 

30 magnitude) a more aggressive value (higher V^^^j) can be computed. 
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Physical amplifiers have finite limits on the rate at which the 
voltage can change. The pulse response calculated above does not include such 
limits. It should be noted that a constraint on the rate of change of voltage in the 
current amplifier can be added to the equation set during numerical solution in the 
5 same manner as all of the other constraints. The slew rate becomes important 
when rapid transients of the input are allowed. The slew rate equation is simple. 
For example, an amplifier chip may lunit at 1.25 volts per msec. The voltage 
constraint used in calculating the input is the minimum of either the voltage 
generated from the Voltage Limit Constraint above or 
10 (old_voltage+(A(I)-A(I-l))* Volt_Rate_Limit) with the correct unit conversions. 

3.2 Rescaling of Current Inputs 

The current conomands generated in the previous section are 
completely scaleable. If for example, a 5 track seek has been created by 
15 generatmg an optimal ctirrent conunand and subsequently shaping it, a 4 track, 
shaped move can be made by simply rescalmg the input by 4/5. The move time 
will be the same as that of a 4 track seek, however, the computation is simple. 
The point is that shaped moves can be scaled in amplitude and they will continue 
to reduce vibrations and acoustic noise. 

20 

4.0 Partial Fraction Expansion Equation System Model 

Once a parameter in the system, such as the fundamental limiting 

parameter, has been identified, an input to the system in terms of that parameter 

can be generated. In order to generate an input which results in reduced system 
25 vibrations, it is first necessary to generate a set of simple equations that defme 

the system. The first step in generating such system equations is to write the 

transfer function for a simple model of the system. 

Figure 1 1 shows a third order model for a disk drive system such 

as that shown in Figures lOA and lOB. For system 80 shown in Figure 11, a 
30 voltage V is input so as to control movement of head 9 to a particular track on 
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a recording medium. To this end, blocks in system 80 convert voltage V into a 
current i, convert current i into torque i, convert torque i into frequency co, 
convert frequency co into velocity v, convert velocity v from meters-per-second 
into tracks-per-second, and convert tracks-per-second into a location of a 
particular track on a recording medium such as a magnetic disk. 

To generate a set of equations for system 80, the transfer function 
H(s) therefor must &st be determined. The transfer function H(s) for system 80 
is as follows: 



K, Ik 

His) = ; , (6) 



where is the motor torque constant, K^is the back EMF constant, 1 is the 
10 actuator arm radius, k corresponds to track density on a recording medium, 1^ is 
the inductance of a voice coil in the motor, 1^ relates to inertia of the motor's 
rotor, R.„ corresponds to an amount of viscous friction in the motor, and 
corresponds to an amount of electrical resistance in the system. 

Next, H(s) is wrinen as an equivalent transfer function in terms of 
15 three parameters, namely A, a and b. as follows: 

where A corresponds to a scale factor and a and b correspond to inverse time 
constants. Specifically, in terms of the parameters in equation 6, A is defined as 
follows: 



• (8) 



em to 



Variables a and b correspond to two solutions in s, namely s, and s 2, to the 
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following quadratic characteristic equation from equation 6: 

IJ + UR^ + R.JJs + RR^ + K,K. = 0. (9) 

em ^ e M e m' em to 

After writing H(s) in terms of the equivalent transfer function 
Heq(s), a modal analysis or a partial fraction expansion is then performed on 
Hcp(s) so that each pole of the system is in a separate term. A partial fraction 
expansion for H«,(s) is follows: 

H(s) = ^* ^ - ^ . (10) 

^ s (-ft + a)is + a) {-a + b)is + b) 

Each term of the expression for H^{s) in equation 10 can then be written as 
follows: 



s 



Ab 



FJs) = ^ (11) 

2^ ' (-b + a){s + a) 

f ^s) - . 

^ (-a + b)[s + b) 

From equations 11, a time domain response of each term to a unit step (e.g., to 
a unit voltage) input can be determined. This time domain response is as follows: 



F,(r) = At 



m - 



Ab 



Ab 



a {-b + a) e'-a (-b + a) 

a (-b + a) 

Aa 



(12) 



Aa 



b {-b + a) e'"b (-b + a) 
Aa . „ 



b i-a + 6) 



(1 -e-»'). 
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System response to a unit impulse can also be obtained by differentiating the 
above "step response** equations with respect to time. 

To solve equations 12 for a system input, values must be 
determined for F,{t), FjO) and F^t). Examination of the expressions for F (t), 
FiCt) and Fjii) reveals that the expressions for F2(t) and F3(t) vanish as "t" 
approaches infinity in a case that the input to the system is zero, which is the case 
at the end time of a move (hereinafter referred to as time "te^d")- Since the input 
is zero at time t^nd. responses of Fjft) and FjCt) decay after time tend- Accordingly, 
in order to solve for a system input that will move the system to the correct 
location in a rest state, FjCt) and ¥^{1) should be set to zero at time t^d. 
Furthermore, since the input is zero after time t^nd, the value of F,(t) stops 
changing at time l^. The above two conditions are equivalent to zero velocity 
and zero acceleration (or current) at the end of the move. 

Thus, in the equations for Fi(t), F2(t) and F3(t), in a case of zero 
velocity and acceleration at the end of a move, only the term involving Fi(t) 
contributes to the final position. An expression can then be written including the 
response of each mode, namely Fi(t), F2(t) and F3(t), at a time corresponding to 
an end of the move, namely x^, for a voltage Vq input at to and havmg transitions 
at ti and tj. To this end, XF/t^ set forth below defines the contribution of Fi(t) 
to the position at time t^nd, ^^^li^nd) set forth below defines the contribution of 
F2(t) to the position at time t^ndi and XF3(t(.„d) set forth below defines the 
contribution of F3(t) to the position at time t^^d 

= ^o^2(^«/ - V - 2 - * 2F„F,a„, - t,) (13) 

Thus, the total movement distance at time t^^ is expressed as 
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XF.iU + XF,{tJi . XF,{t^ . (14) 

As noted above, the expressions for F2(t) and F^t) vanish as **t" approaches 
infinity in a case that the input is zero, i.e., at the end of a move. Thus, in 
equations 13, only the term involving Fi(t) contributes to the final position as time 
approaches infinity. This results in the following three constraint equations that 
5 must be satisfied in order for a third order system to reach its desired final 
position in substantially minimal time, with substantially zero velocity and 
acceleration 

^\i^end) = Firialpos 

^2i^eJ = 0 (15) 

where Finalpos denotes the end (or "final") position of the system. These 
equations correspond to rigid mode motion constraints since they achieve a final 

10 position in the third order model with substantially zero velocity and acceleration. 

A digital model of the response of a disk drive system to a unit 
impulse may be derived based on the foregoing. Figure 24 shows a digital 
voltage input waveform. The waveform of Figure 24 can be conceptualized as 
a series of pulses as shown in Figure 25. Each of these pulses can be represented 

15 as a pair of impulses, as shown in Figure 26. A first digital constraint is formed 

by realizing that the system^s final position, Finalpos, is proportional to the area 
enclosed by the voltage waveform, i.e., the integral of the voltage waveform. 
The area of the voltage waveform can be determined by inspection of the 
waveform shown in Figure 24. Specifically, the area is determined as follows 

A' 

Finalpos = X) ^/^ > ^^^^ 

20 
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where "VjAt" corresponds to an area of a pulse "i". Negative values of V ^ 
subtract from the overall sum. 

The remaining two digital constraints can be formed using the 
impulse summation form of the voltage waveform shown in Figure 26. 
Specifically, as noted above, the impulse response of the system can be 
determined by taking the derivatives of step response equations. For FjCt) and 
F3(t) above, the impulse response equations are as follows: 



dt 



= f;^ 



(17) 



In accordance with constraint equations 12, the response of the series of impulses 
of the voltage sunmiation form of the voltage waveform, evaluated at time tend, 
is set to zero in order to form the second and third digital constraints, as follows: 



/=1 

N 



1=1 ' b - a 



(18) 



/=i ' a " b 

Thus, partial fraction expansion equations corresponding to a digital model of the 
response of a disk drive system to a unit impulse comprise: 
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Finalpos - V^AAt 



1=1 





(19) 



where Finalpos, A, a, b and t^nd are as previously defined, Vj are voltage inputs 
to the system, are the times at which Vi are input, and At is a time interval at 
which Vi are input. 



in substantially zero velocity and acceleration at the end of a move. These 
equations can be solved in a number of ways. For example, the equations can 
be solved using a linear programming method, either on their own or in 
conjunction with a variety of other constramts. Examples of such constraints 
include vibration constraints and physical constraints such as power consumption 
or heat buildup, etc. Also, the continuous form of the equations can be solved 
for voltage switch tunes, as described in more detail below. 

5.0 Generating "Good Moves'* 

The techniques set forth in sections 3 and 4 above determine inputs 
to a system that utilize substantially the full capabilities of the system acmators 
and that reduce system vibrations (both mechanical and acoustic). Several 
techniques can be used, either alone or in conjunction with the techniques 
presented above, to generate (or enhance) system inputs to reduce system 
vibrations even further. 



defined as a movement of a system component along a trajectory from an initial 
position to a nominal final position, such that, at the time the component reaches 
its nominal final position, the component is vibrating at or below an acceptable 



Equations 19 can be used to solve for system inputs Vj which result 



For the purposes of the subject application, a "good move" is 
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level, possibly subject to one or more predetermined constraints. In the context 
of computer disk drives, this means that, as soon as the drive's head reaches its 
nominal final position, the head is ready to perform a read/ write operation to a 
recording medium. A good move is also one in which acoustic noise excited by 
the move is at or below an acceptable threshold. In contrast, a "bad move" is, 
for example, a move in which the component is vibrating at an unacceptable level 
when the component reaches its nominal final position. In the context of 
computer disk drives, this means that at the time the drive's head reaches its 
nominal final position, the head is vibrating at a level which is too high to 
perform a read/write operation accurately. A bad move is also one in which the 
acoustic noise excited by the move is above an acceptable threshold. 

A powerful result of the approach outlined herein is that a rigid 
body conunand may be calculated by one of a number of techniques including 
equations 19 or the current mode technique of section 3.1. Next, other dynamics, 
such as flexible modes, acoustic modes or other dynamic effects can be 
subsequently corrected by the various methods set forth below. For disk drives, 
the process comprises, e.g., determining, over a population of drives and for each 
move length (or range of move lengths), which frequencies can be heard and 
which effect settling time. Then, moves may be generated that limit the 
amplitude of those fi-equencies or some subset thereof. Frequencies that only 
appear in a subset of the drives may also be included or tumed-on selectively. 
The effects of the selected frequencies can be reduced by applying any of the 
techniques set forth below. 

5.1 Vibration Limiting and Sensitivitv Constraints 

This technique for generatmg a good move incorporates physical 
constraints into a system model. One such model comprises partial fraction 
expansion equations 19 derived above. To generate a good move in accordance 
with this technique, the first step is to choose one or more physical constraints. 
Examples of physical constraints include constraints that limit vibration (both 
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mechanical and acoustic) and constraints for insensitivity to parameter variations. 

Constraints for insensitivity to parameter variations can be added 
to the system equations in several ways. For example, derivatives of the partial 
fraction expansion equations can be taken with respect to a particular parameter 

5 and then added to the constraint equations. This technique is set forth in detail 
in U.S. Patent No. 5,638,267, issued to Singhose et al. on June 10, 1997, the 
contents of which are hereby incorporated by reference into the subject 
application as if set forth herein in full, and in U.S. Patent No. 4,916,635, issued 
to Singer et al. on April 10, 1990, the contents of which are also hereby 

10 incorporated by reference into the subject application as if set forth herein in fiill. 
Additionally, a sampling technique can be used to set a series of constraints for 
different values of a given parameter so as to limit the solution for a variation in 
that parameter. For example, using this technique, 10% insensitivity at 1000 Hz 
can be achieved by setting constraints at 900 Hz, 920 Hz, etc. up to at least 1100 

15 Hz. One variation of this technique is set forth in U.S. Patent No. 5,638,267. 

In accordance with the invention, constraints can be chosen so as 
to reduce system vibrations by substantially 100% or to limit reduction in system 
vibrations to a percentage which is less than 100%, say, between 50% to 70%. 
Which constraints are used to reduce system vibrations by tuch amounts can be 

20 determined experimentally or by running a computer simulation for the system. 
Choosing constraints that reduce system vibrations by less than 100% (e.g., a 
range of 50% to 70%) is particularly useful m cases where excessive vibrations 
do not occur in the system and in cases where large reductions in vibration?, 
including noise-inducing vibrations, are not required. In these cases, since the 

25 constraints are more relaxed, the resulting system input is shorter in duration. 
As a result, the system moves faster than it would otherwise. In general, there 
is a trade-off between move time and the use of stringent constraints. That is, 
while the use of more stringent constraints may further reduce vibrations, move 
time is slowed as a result of their use. 

30 Vibration constraint equations may also need to be scaled in order 
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to achieve desired vibration reductions. Specifically, in terms of the foregoing 
model of a disk drive system, equations for vibradon amplitude are relative. This 
is because the rigid mode equations do not take into account flexible body 
equations (described below in section 6.3). If a full set of equations (i.e., both 
rigid and flexible body equations) were derived for the system, and then a partial 
fraction expansion applied thereto, the coefficients of the vibrational terms would 
scale to the problem. However, since the equations were derived separately and 
later combined using the techniques described herein, the relative scaling is not 
known for the equations given above, and the problem must be scaled manually. 
In the case of zero vibrations, no scaling is necessary. For any other case, one 
of several scaling methods may be used. For example, an arbitrary gain may be 
selected and a resultmg sequence applied to the system. The resulting vibrations 
can then be measured and the vibration constraint equations multiplied by a ratio 
of a predicted vibrational level divided by an actual measured vibrational level. 

By way of example, consider the fourth order transfer function 
shown below in equation 20. The transfer function corresponds to the position 
response of a simple mass-spring-mass system to a force input. 

His) = 4 * ~ T (20) 



The first term of the transfer function, namely 

C 



2 - (21) 



corresponds to the rigid body motion, while the second term, namely 



(22) 



+ 2zb>S + 00^ 
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corresponds to a vibratory mode with undamped natural frequency "co" and 
damping ratio "z". When transfer function 20 is expanded in a panial fraction 
form, the following three terms result: 



//(,) = Z2££ . £ . C(-a) . 4z^o) . 2zs) ^23) 



The third term from equation 23, namely, 
C(-o) ^ 4z^o) + 2zs) 



(24) 



corresponds to the transfer function of the vibratory mode, and should be 
compared to the second order vibratory term included in the original transfer 
function. Using inverse Laplace transformation, the time domain impulse 
response of the vibratory mode is determined as 



-Cg'-'^^siDh(o)v^-l + t) ^ 



2Cz^e'^'siDh(a)^"l + z^ t) ^ 



(25) 



2Czg-^"^cosh(o)v/"l ^ t) 



0) 



Using suitable trigonometric identities, the time domain response of the 
vibratory mode in equations 25 can be simplified as follows: 

-Ce --"'sinCw^O ICz'e --'^^'sinCco^ 0 ICze -*'"'cos(co^O 



+ 



+ — —, (26) 



or, after collecting terms, to: 
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Ce * [!l i ^ + ^] , (27) 



where "co<," is the damped natural frequency 

CO, = a)/r^. (28) 

The impulse response expression derived above can be used direcrty to form a 
quantified constraint equation; i.e. the magnitude of the constraint equation can 
be set to be less than or equal to an actual vibration level in position. 

It should be noted that the general form of equation 27 is a phase 
shifted sine wave of frequency coj. The sine wave is further subject to a decaying 
exponential envelope with time constant 



— . (29) 



Equation 27 can be rewritten again, including the phase shift ^ and scale factor 
C,, as follows: 

Cie--^sin(a)^r + 4>) . (30) 

This form of equation 27 should be compared with the expression 
for the impulse response of an isolated second order system (as defined by the 
second term of the original transfer function H(s) 



_ — — - ) (31) 

s + 2zu)s + <i)" 



-36- 



wo 99/45535 



PCT/US99/04910 



shown below: 



— *e'=^'sm(o).r + 4)). (32) 



Since we are only concerned with the magnitude of the residual vibration, phase 
shift $ is of no hnportance, and the expression for the isolated second order 
system can be used in the formulation of quantified constraints if the scale factor 
co/coj is replaced by Ci- Parameters Ci, zo) and cOd can be determined 
experimentally by (i) subjecting the actual system to an impulse input, and from 
experimental data, (ii) estimating the period of vibration and computing cj^ as 
follows 



co^ = — , (33) 

T 



10 where i is the period, (iU) estimating the decay envelope to determine zco, and 
(iv) estimating actual amplitude to determine Cj. 

After constraints have been chosen based on the foregoing criteria, 
a trajectory is digitized and is then solved for each trajectory command value in 
time subject to the chosen constraints. Digitization of the trajectory results in 

15 simpler, and often linear, equations that can be solved readily for the complete 
system trajectory. The invention, however, is not limited to solving the system 
equations in this manner. That is, any known technique for solving the system 
equations can be employed in the context of the present invention. The 
digitization rate used is a factor that contributes to the successful implementation 

20 of the invention. This issue is discussed in section 7.1 below. 
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5.2 Input Shaping™ 

This technique for generating good moves uses Input Shaping'^^ and 
its substantial equivalents. One example of Input Shaping^^^ comprises 
convolving a function with a system input in order to reduce unwanted system 
vibrations. Examples of system inputs include current and voltage. Specific 
Input Shaping™ methods are set forth in detail in U.S. Patent Nos. 5,638,267 and 
4,916,635, which were incorporated by reference above. 

In accordance with this technique, an input to a system, such as a 
voltage command or a velocity, is shaped in order to produce a desired output. 
For example, the input may be Input Shaped'^'^^ by convolving the input with an 
impulse sequence which reduces unwanted vibrations in the system, including its 
flexible modes. Alternatively, the input may be shaped using any one or more 
of the filtering techniques presented herein, or using any other known conmiand 
profiling techniques. 

5.3 Filtering 

In addition to Input Shaping^*^, as noted above, any type of 
conventional shaping method can be used to generate an input that results in a 
good move. For example, it is also possible to generate a good move by filtering 
a trajectory with infinite impulse response ("IIR") or finite impulse response 
("FIR") filters, including filters where there are zeros near poles of the system, 
such as those described below. 

5.4 Least Squares Solution to Frequency Sampling 

A conventional least squares approach can also be used to generate 
an input that results in a good move for a computer disk drive system. An 
example of such an approach is provided in U.S. Patent No. 5.638,267. This 
approach may also be used to limit high mode excitations. 
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5.5 Selective Shaping'^'^ 

Input Shaping""^, filtering, and the like can increase seek time in 
some cases. However, the length of an Input Shaper™ is directly related to the 
period of the vibration being reduced. At low frequencies, i.e., long periods, the 
increase in seek time caused by an Input Sh^er'^^ can be a substantial percentage 
of the total seek time. This is of particular significance in relation to computer 
disk drives. That is, in computer disk drives, relatively low frequency vibrations 
can occur due to compliance of the drive's mounting in its enclosure. For 
example, a 200 Hz mounting resonance will have a period of 5 ms. In this case, 
a three-pulse Input Shaper™ will also have a period of 5 ms. Thus, if the 
computer disk drive has an average seek time of 10 ms, the drive's average seek 
time will increase by 50% if the Input Shaper'^'^ is applied to all moves. 

In view of the foregoing, it is advantageous to apply Input 
Shaping™ selectively. This is possible because not all moves result in unwanted 
vibrations. That is, some moves are inherently good, meaning that they produce 
relatively few vibrations, and some moves are inherently bad, meaning that they 
produce excessive vibrations. By evaluating constraint equations for a given 
move, whether that candidate move is good or bad can be predicted in advance. 
That is, a move can be evaluated to determine whether thcmove excites greater 
than a predetermined level of vibrations in the system. Input Shaping™ can then 
only be applied to moves that excite greater than the predetermined level of 
vibrations. Thus, only those moves exciting greater than the predetermined level 
of vibrations, e.g., the bottom 30% of the moves, need be Shaped ™. 

By virtue of the foregoing, the present invention provides the 
advantages flowing from using Input Shaping'^^ with less of an increase in seek 
time than would otherwise be achievable. Moreover, in accordance with the 
invention, a level of Input Shaping™, or- the numbers and types of constraints, 
may be altered depending upon how good or how bad a move is, i.e., depending 
upon the level of vibrations that the move excites. For example, if a move only 
requires minimal Input Shaping^'^', lengths of impulse sequences used to Input 
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Shape^^ the move may be reduced. 

In the context of computer disk drives, for a given frequency of 
interest, a measure of the goodness of a move can be evaluated by examining a 
superposition of components of an input at the frequency of interest. An example 
of such a superposition is shown in the graph of Figure 12, in which the 
frequency of interest is 200 Hz, and in which a good move/bad move threshold 
is labeled 81. By examining threshold 81, it is evident that movement distances 
within a range around 1800 tracks and 7000 tracks are inherently good moves. 
On the other hand, movement distances requiring vibration reduction are 
identified in Figure 12 by cross-hatching. For these movement distances, the 
amount of vibration reduction required is determined by taking a ratio of a 
move's goodness level to the good/bad threshold. The magnitude of the Input 
Shaper*^ required for each movement distance is determined based on this ratio. 

Specifically, the magnitude of an Input Shaper™ response can be 
equated to a desired reduction in vibrations ("mag") to derive an expression for 
a modified Input Shaper'"^. The expression for the period, t , of a modified two- 
pulse Shaper™ is as follows 

^ ^ arccos(2mag^ - 1) ^^^^ 

27C/ 

It is noted that the above expression is closed form and, therefore, is suitable for 
on-line, real-time implementation. For unity magnitude Shapers^"^ (mag=l) in 
which no vibration reduction occurs the foregoing expression yields 

T = 0, (35) 

and the two pulses are co- incident at time t=0. For zero magnitude Shapers"^^* 
(mag=0) in which 100% vibration reduction occurs, the foregoing expression 
yields 
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-i. (36) 



or '/2 period. At intermediate values of mag, the expression yields intermediaie 
values of the Shaper""^" period x. For example, when mag=0.5, 



T = i- . (37) 
3/ 



At this point, it is noted that although two and three impulse Input 
Shapers'^'^ are described above, the invention can use Input Shapers '^'fiaving 
5 more than three impulses as well as to shapers designed using other techniques. 
It is also noted that although the invention has been described with respect to 
Input Shaping™, any shaping technique may be used. For example, moves that 
are considered bad moves could be selectively shaped using an FIR filter, an IIR 
filter, any shaper, or even just extended in some arbitrary way to achieve the 
10 goals of the present invention. In addition, constraints, such as vibration limiting 
constraints or the like, may be selectively added to the system equations based on 
whether a move has been identified as good or bad, and the resulting set of 
equations can be solved for a complete move trajectory. The moves also may be 
solved off-line and tabulated to achieve the same benefit, 

15 

5.6 Voltage Control Bv Forcing Switch Times 

Disk drives operate in a current command mode, meaning that their 
operation is controlled by current commands. Since saturation is in voltage, 
conventionally there is no way to determine if a particular current command will 
20 saturate the system (notwithstanding section 3.1 above). Accordingly, 
commanding the system in current can fail to produce a desired movement. The 
invention addresses this problem by commanding a current that is guaranteed to 
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saturate a voltage input to the system. In this way, the inventioQ is able to 
command voltage in either direction. That is, it is possible to provide current 
commands to the system, but actually to be commanding voltages switches. 

In more detail, many disk drive systems utilize a closed-loop 
current control system to control current through a voice coil motor. A closed- 
loop controller in such a system makes it possible to specify a current command 
that, when applied to the motor, results in a constant torque. Specifically, the 
closed-loop controller adjusts the voltage through the motor to maintain a desired 
current level. However, there are limitations to such a controller. Specifically, 
the maximum voltage available to the controller is limited by the power supply 
used in the disk drive. Thus, if the voltage required to achieve a particular 
current is above the power supply's maximum voltage, the controller will saturate 
at the maximum voltage and the desired current will not be achieved. 

While the foregoing may be perceived as a drawback, it can 
actually be an advantage, since the saturation can be used to operate the current 
controller in an alternative mode m which power supply voltage is simply applied 
to the motor's terminals. That is, the controller is operated as an open-loop 
voltage switching system, instead of as a closed-loop current control system. As 
a result, it is possible to create an alternative formulation for control strategies 
and motion profiles. This technique of using current control to command voltage 
has the additional benefit that voltage conmiands can be used during disk drive 
seek; however, a normal current conmiand mode is still available for tracking 
without changing any system hardware. 

Use of the open-loop voltage control technique described above 
may require a series of precisely-timed commands to achieve vibration 
suppression effects. Section 7.1 below presents a variety of techniques for 
implementing such precisely-timed commands. If the commands are not in the 
fundamental limiting parameter of the system, then section 3.1 may also be used 
to generate such commands. If the commands are in the fundamental limiting 
parameter, then the techniques of sections 4 and 5 can be used. 
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In this regard, the fundamental limiting parameter is often not used 
to conunand the system due to hardware constraints (e.g., the tracking controller 
in a disk drive system is designed for current, not voltage, commands). In the 
case of a disk drive, by modifying the drive's amplifier, it is possible to use 
voltage conunands during one mode of operation, e.g., drive seek operations, and 
to switch to current commands during another mode of operation, e.g., drive 
tracking operations. Thus, for example, the techniques of section 3.1 can be used 
during a tracking operation and the techniques of sections 4 and 5 can be used 
during a seek operation. 

5.7 Other Techniques 

In addition to the foregoing techniques, other techniques may also 
be used to generate (or enhance) an input to produce a good move. For example, 
it is possible to adjust any of the above techniques using real-time identification 
of any system parameters. Specifically, a system parameter, such as temperature 
or resistance, may be identified, and that parameter addressed when generating 
a good move. Additionally, symmetric and muUi-step (i.e., Posicast) inputs can 
be used to generate a good move. As set forth in section 11 below, transient 
portions of the input or feedforward signal at the beginning, middle and/or end 
of a move may also be varied in order to reduce vibrations resulting from the 
move. Additional constraints may also be added to the system in order to control 
vibrations. For example, these constraints can be symmetric and can be a 
function of time or position and/or can relate to thermal limits, current limits, 
duty cycle, and/or disk drive ZV constraints. Moreover, it is noted that since 
each transition of an input corresponds to a complete move in and of itself, any 
of the above techniques for generatmg a good move set forth above can be 
performed in real-time, and the Shaping™ and filtering methods described above 
can be performed only on the rapid transient portions of the move. See section 
11 below. 
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The Shapers™ and filters described above are preferably stored in 
a memory as computer-executable process steps, and executed by a processor in 
the system. However, the invention can also be implemented in a fuzzy logic 
controller by embedding either an Input Shaper™ or direct solutions to systems 
5 of constraint equations in a fuzzy logic block. Sunilarly, the invention can be 
implemented in a neural network controller that can be constructed by imbedding 
either an Input Shaper™ or direct solutions to systems of constraint equations into 
a neural network. Both the fuzzy controller and the neural network PV controller 
can be designed to handle non-linear system parameter variations which enter into 
10 equations for a system. 

6.0 Controlling Trajectory as a Function of Time 

This aspect of the invention is directed to a control technique for 
reducing unwanted vibrations in a dynamic system by controlling a trajectory 

15 which is specified as a function of time. More specifically, in this aspect of the 
invention, an input that corresponds to a good move is used as feedforward 
trajectory for the system. Any of the above techniques for generatmg inputs 
corresponding to a good move may be used to modify inputs in the following 
control system. However, for the sake of brevity, only Input Shaping"^ will be 

20 described in detail. Moreover, although the invention can be used to control any 
dynamic system, a description will be provided only for a computer disk drive 
system, again for the sake of brevity. 

6.1 Optimization 

25 One approach to generating an input that will result in a fast 

response in a computer disk drive is to perform an optimization using partial 
fraction expansion equations 19. To perform the optimization, these equations 
are augmented by including a set of constraints to address vibrations in various 
modes of the system. In addition, constraints may be added for vibrations m 

30 modes that are near a predicted vibrational mode and for real world conditions, 
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such as voltage limits. A digital fonnulation of the problem can then be modeled 
so as to make all of the constraints linear, thereby simplifying their solution. The 
digitization rate may be selected using the techniques presented in section 7.1. 

The invention is not limited to methods for performing optimization 
described herein. Rather, any optimization technique may be used. . For example, 
other methods for implementing optimization which can be used in conjunction 
with the present invention include off-line computation and use of efficient linear 
solving codes on the disk drive controller. It is also possible to generate a series 
of solutions off-line and to fit curves to the solutions. Appendix C shows a set 
of linear equations written in GAMS that can be solved for ZVD constraints on 
a problem in which voltage is the commanded input. 

6.2 A Related Approach 

The partial fraction expansion model of a disk, drive system shown 
above can also be used in conjunction with Input Shaping'^^ in order to generate 
system inputs. For example, in an article entitled "Fast Servo Bang-Bang Seek 
Contror by Hai T. Ho, ®1997, IEEE, the contents of which are hereby 
incorporated by reference into the subject application, it was demonstrated that 
near the end of a move, the instantaneous state of a system can be measured. 
Inputs to the system can then be calculated so as to bring the system to its final 
location with limited vibrations. The flaw in this technique, as presented in the 
article, is that exact solutions to the disk drive system equations could not be 
obtained because a second order system model was assumed. Since a disk drive 
system is more accurately modeled as a third order system, this technique could 
result m inputs which did not sufficiently limit vibrations or which resulted in the 
system reaching a nominal final position without conung to rest. By use of the 
partial fraction expansion model presented in section 4, a more precise solution 
to the system equations can be obtained, thereby resulting in relatively fewer 
vibrations when the system reaches its final location. Moreover, insensitivity 
constraints can be added to the foregoing equations by a variety of techniques, 
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e.g., by sampling (as discussed in section 6.7), by taking derivatives (as discussed 
in U.S. Patent No. 4,916,635), or any other approach, thereby leading to robust 
solutions of the equations and still fewer resulting vibrations, 

6.3 Rigid and Flexible Mode Problem 

Another approach which uses Input Shaping™ in conjunction with 
the partial fraction expansion model of section 4 is to separate the problem into 
a rigid mode problem and a flexible mode problem. Examination of modal terms 
associated with a disk drive reveals terms similar to a third order system, e.g., 
F,(t), FjCt) and F3(t) from equations 12, plus additional terms associated with 
oscillating or vibrational modes, e.g., F4(t), FgCt), etc. (not shown). The present 
invention separates these terms into two sets. The first set, which defines the 
rigid mode of the system, includes the third order system modal terms, i.e., Fi(t), 
F2(t) and F3(t). The second set, which defines the flexible modes of the system, 
includes the other modal terras, i.e., F4(t), F5(t), etc. The partial fraction 
expansion equations in section 4 can be solved for inputs which drive the system 
so as to satisfy rigid body motion constraints F,(t), F2(t) and F^ii) and thereby 
achieve a final position with substantially zero velocity and acceleration. 
However, vibrations resulting from the other modes, namely F4(t), F5(t), etc., can 
still be problematic. 

More specifically, the partial fraction expansion equations cannot 
be solved for an input which reduces vibrations in the flexible modes without 
adding additional constraints for each additional mode. However, adding 
additional constraints increases the complexity of obtaining the solution. 
Additionally, if curve fit or empirical solutions are to be obtained, the solution 
is increasingly more complicated for each additional constraint that is imposed. 
Therefore, in order to reduce vibrations in the flexible modes without 
substantially increasing the complexity of obtaining a solution, die present 
invention first solves the rigid mode equations in order to determine system 
inputs. Next, the determined system inputs are shaped using an Input Shaper"™ 
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designed for the flexible modes of the system. For example, the Input Shaper'^M 
may be designed based on frequency and damping ratios of F4(t), etc. above. 
The partitioning of Fj(t), F2(t) and F^{i) into the rigid mode and 
F4(t), F5(t), etc. into the flexible mode is merely an example of one way to 
5 partition the system's modeled modes. The invention is not limited to such 
definitions. In fact, the partitioning of modes can vary depending upon the 
system under consideration as well as various other factors. For example, it is 
possible to derive a second order system input based on a system model that 
includes only Fi(t) and F2(t) in the rigid mode and that includes FjCt), F4(t), etc. 
10 in the modes that are addressed using Input Shaping™. Such a formulation would 
be of particular advantage in the case where F3(t) corresponds to a system's 
electrical time constant, since Input Shaping^^^ could be used to correct variations 
in the system's electrical time constant. 

15 6.4 Center Of Mass 

In the rigid mode of a system, problems may arise because various 
components of the system are moving relative to one another. For example, in 
the case of a computer disk drive, vibrations may be induced in the system 
because the motor's rotor is moving out of phase with an endpoint, such as the 

20 read/write head. This will cause problans during feedback if a rigid body model 
is assumed. That is, if die endpoint alone is used as feedback, the deflection that 
is observed will be fed back and cause corrections that are not required and thus 
may be detrunental to system performance. To reduce such errors, the present 
invention estimates the center of mass of the system and uses the estimated center 

25 of mass, rather than an endpoint, as the system's feedback. One way to estunate 
the center of mass of the system is to calculate a state estimator using modern 
control techniques. A wide range of correction factors, such as a quasi-static 
correction factor similar to that described below, can be applied to improve the 
estimate of the center of mass coordinate. 
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The center of mass coordinate estimated above will not reflect the 
vibrations of the rotor nor those of the endpoint. Consequently, the rigid body 
behavior of the system can be decoupled from most all vibrations (including 
acoustics) in order to allow for more accurate feedback control than has 
heretofore been possible. 

6.5 Quasi-Static Correction 

On the other hand, if endpouit feedback, rather than center-of-mass 
feedback, is used in the system, vibration-reducing inputs to the system, such as 
those generated in accordance with the techniques set forth herein, can cause 
unusual side effects. That is, the component at the endpoint (e.g. , the head in the 
case of a disk drive) does not vibrate. Rather, the component deflects during the 
move. The deflection is in one duection during acceleration phases and in 
another direction during deceleration phases. The feedback signal to the 
controller is measured based on this deflection. As a result, if the system 
controller uses this feedback signal to compare to a trajectory, the controller will 
detect an error even if none exists. To correct this problem, a quasi-static 
correction factor may be added to, or subtracted from, the feedback signal 
depending on the amplitude and the direction of acceleration. 

6.6 Closed-loop Control 

By adopting a control strategy in which a trajectory is specified as 
a function of time, a closed-loop system formulation can be made. As shown, for 
exanq)le, in Figure 11, feedback loops may be closed around any combination of 
system states. In the case of a disk drive, these states include voltage, current, 
acceleration, velocity, and position. An mput which conesponds to a good move 
can then be generated based on the closed-loop formulation. This "good" input, 
which, in the case of a computer disk drive, may be a voltage input, is used to 
generate feedforward signals to some or all of the closed-loops. 
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Feedforward trajectories for each of the parameters to be used as 
a feedforward term are derived from a trajectory selected as a "good" move. 
The technique for generating feedforward trajectories may be based on either a 
simple simulation of the system or the partial fraction expansion model set forth 
in section 4. When generating a good input, it may be desirable to compute a 
move assuming that less voltage (or whatever physical variable is limited by 
reality) is available. This increases the chances that the system will stay on a 
trajectory corresponding to the good move. Another approach to generating a 
good input is to assume full voltage (or whatever physical variable is limited by 
reality) and to assume safe values for the physical system parameters used in the 
model. In this case, full vohage would only be used if the parameters were the 
worst case possible. At all other nominal cases, a lower voltage would be used. 

In the closed-loop control model, the feedforward terms for each 
loop could be replaced by only one feedforward (or input) to the outer loop, 
where this one feedforward term is computed as a combination of each of die 
other feedforward terms. Also, feedforward inputs for the system could be scaled 
differently for each system parameter. For example, different gains could be 
used for each of velocity, position, acceleration, current, and voltage feedforward 
terms. The selection of gains adds zeros to each feedforward term. Since Input 
Shaping*"^ is a method of placing zeros in certain locations using a constraint 
based technique, it is understood that the gains selected for the feedforward terms 
can be chosen to add zeros which serve to augment Input Shaping™. This 
technique would lower the requirements on any Input Shaper™ used, or enable 
fewer constraints to be used in the derivation of complete trajectories. 

Shown in Figure 27 is a block diagram of a simple system, namely 
system 82, that will be used to demonstrate how a collection of feedforward 
signals can be used to effect zeros in a transfer function. In system 82, the 
variable to be controlled is the output position x, there is a feedforward input FF, 
and there are three command inputs, namely position input x^, velocity input v^, 
and acceleration ii^ut a^. There is also full state feedback of position x, velocity 
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V, and acceleration a. In the ideal case, the feedforward input FF will cause the 
system to follow a desired trajectory (x, v, a) through state space. If the actual 
trajectory does not match the desired trajectory, then error terms (x^ - x), (v^ - v) 
and (a^ - a) will be multiplied by gains K,, and to generate corrective 
inputs. If the inputs v^,, a^ and FF are all set to zero, die system will still "track" 
the remaining conunand input x^, but likely with larger position error. 

The command inputs, x^, and a^,, and the feedforward input FF 
are produced by another process and are simply inputs to the system. However, 
the command inputs x^, v^ and a^ are mathematically related as derivatives. Their 
relationship is indicated in the block diagram shown in Figure 28 by applying the 
Laplace differentiation operator "s" repeatedly to the x^ input signal. It is well 
known that it is physically impossible to implement the zero corresponding to the 
"s" operator without a corresponding pole. Since the input signals Vj and a^j are 
actually derivatives, the use of the "s" notation in Figure 28 is a notational 
convenience facilitating intuition, not a means of implementation. 

By applying block diagram algebra to the block diagram shown in 
Figure 28, scaled versions of all the inputs (Xd, v^, a^ and FF) can be shown to 
enter the closed-loop system at a single point, as shown in Figure 29. The 
combination of the three parallel inputs (to the nominal closed-loop system) 
shown in Figure 29 can be reduced further using block diagram algebra as shown 
in Figure 30. The numerator of the block shown in Figure 30 is a third order 
zero. For this system, the zeros associated with the combination of feedforward 
terms exactly match (and cancel) the poles of die nominal closed-loop system. 
Thus these zeros serve to cancel the dynamics of the modeled system. The 
location of the zeros can be adjusted by including weighting factors Wy, W. and 
Wf on each of the terms v^, a^ and FF, respectively. Following through with the 
analogous block diagram reduction yields the result shown in Figure 31. 

By allowing additional freedom in the placement of the zeros (by 
adjusting weights W^, and Wf), ±e zeros can be placed to effect a first order 
cancellation of the (un-modeled, but measured) poles of the system, e.g. poles 
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associated with vibration or other unwanted dynamic response. Since FIR input 
shapers can be viewed as "all zero" filters, the zeros associated with the weighted 
feedforward input can be used to replace or augment the action of the zeros of the 
FIR input shaper. Thus, since some zeros are supplied by the weighted 
5 feedforward input, a lesser number of zeros must be supplied by the shaper, and 
consequently, the shaper period becomes shorter. The overall time advantage lies 
in the fact that the zeros associated with the weighted feedforward are 
implemented without a tune cost (i.e., without a penalty). 

10 6.7 Additional Constraints 

It is also possible to add constraints to the partial fraction 
expansion equations of section 4 in order to provide for relatively rapid 
movements of the rigid mode portion of the system. For example, it is possible 
to add constraints relating to insensitivity to variations in system parameters. To 

15 do this, one or more constraint equations in the continuous form of equations 19 
may be differentiated with respect to the system parameter and the resultant set 
to zero. This ensures that the solutions to the equations are insenshive to 
variations in that parameter. A similar result can be obtained by solving the 
equations to which insensitivity is required multiple times .with different values 

20 for the system parameter. For example, the second of the above partial fraction 
expansion equations 19 can be repeatedly solved with different torque constant 
values in order to assure that the solution thereto works well for a range of torque 
constant values. Such an equation must be given some degree of freedom by not 
requiring that the solution be exactly zero for all values of torque constant. 

25 

6.8 Switch Times 

As noted above, the continuous form of the partial fraction 
expansion equations can also be solved for voltage switch times which reduce 
vibrations in the system. A voltage switch time corresponds to a transition time 
30 between full positive and full negative voltage inputs to the system. In solving 
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for these switch times, if no other constraints are included, the system can be 
modeled using the three partial fraction expansion equations 19 shown above in 
section 4. Using the continuous version of these equations, it is possible to solve 
for three transition times, namely t^, to and tend* as set forth below. 
5 Equations 19 are nonlinear in their continuous form. However, 

they can be manipulated to make them easily solvable in a controller with limited 
computational power. Since the first equation is linear, 'the first equation can be 
substituted into the other equations, yielding two equations with two unknowns. 
The solution can then be obtained by performing a simple fixed-point iteration. 

10 Only a few iterations are needed to obtain a solution to machine precision. 
Additionally, the starting guess for the iteration can be obtained from a simple 
curve fit, thus further increasing the speed of the calculation process. It is also 
possible to generate solutions to these equations by defining switch times as a 
function of movement distance and then fitting the data with an empirical formula 

15 to achieve a solution that is sufficient for use in a digital controller. 

To this end, an additional step that is performed in order to use the 
continuous switch time solution in a digital controller is to round-off the switch 
times to a nearest digital time spacing. However, because significant errors can 
be generated by rounding-off the switch times, another approach is often used. 

20 The switch times are rounded upward and the voltage levels in each region are 
adjusted to achieve a solution to die same equations. Alternatively, voltage levels 
in only a portion of the solution are adjusted in order to achieve the same result. 
For example, only a few digital outputs are adjusted to correct for the round-off. 

The pseudo-code shown in Appendbc A illustrates (i) the calculation 

25 of continuous switch times, (ii) a method of rounding the continuous switch times 
to discrete time intervals, and (iii) recalculation of the voltage levels to satisfy the 
original constraint equations including discrete switch times. In the first section 
of the pseudo-code, the physical system parameters {V^^, K^, K^, 1, k, Re, R„„ 
and J„j) are used to compute transfer function parameters (A, a, b). While in 

30 theory this step is not required (a formulation in the original system physical 
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parameters could be used), the resulting equations are much more compact and 
more suitable for real-time implementation. 

In the formulation shown in Appendix A, the accommodation of 
digital switch times requires that the entire voltage level of each of the three steps 
be modified in order to satisfy the constraint equations. Modification of the 
entire step resulted because of the specific formulation in terms of three steps in 
voltage. If the problem is reformulated, e.g. in discrete impulse response form 
corresponding to partial fraction expansion equations 19 (repeated below as 38 for 
convenience), other, more flexible modifications of the voltage wave form are 
possible. 

AT 

Finalpos =53 V^Atd 

0 =X) V^-^ie"'^^-''''^'^^ -e'"^'''^'''^ (38) 
f»i b " a 

/=! * a " b 



The pseudo-code shown in Appendix A can be used to compute the 
digital switch times (rounded from continuous), and the impulse response form 
of the constraint equations used only for the final adjustment of voltage levels. 
A voltage variable for the reformulated voltage adjustment problem can be 
formed by selecting an individual Vj, or alternatively selecting a dependent set of 
Vf (e.g. Vi=V2=V3=Vfi„J. The three constramt equations form a linear set of 
three equations in the three variables and can be solved using any method suitable 
for linear equations. If more than three variables are chosen, then the solution 
will not be uniquely determined and opportunity exists for additional criteria to 
be used to select from among the possible solutions. 
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6.9 Shaped^^' Switching 

According to the well-known bang-bang principle, time-optimal 
system commands for a voltage-limited system like a computer disk drive 
comprise a series of full-positive and full-negative voltage commands. In such 
a case, the voltage levels are known, but their optimal switch times are not. 
Because optimal switch times for the system are unknown, it is likely that 
unwanted vibrations will be excited when the voltages commands are input to the 
system. Several techniques for reducing unwanted system vibrations can then be 
performed. For example, it is possible to determine a switching sequence 
analytically subject to various constraints, such as vibration constraints, 
insensitivity constraints, or any number of unwanted dynamic or physical system 
constraints. Additional techniques that can be used include determining switch 
times using any technique described herein and adjusting the switch times in real- 
time to take system movements into consideration. 

Another effective technique for commanding a disk drive to move 
is 10 combine the benefits of a bang-bang time-optimal conmiand with the benefits 
of Input Shaping™. To this end, switch times for an input to a rigid system are 
determined analytically, and then an Input Shaper™ is applied to the input in 
order to shape it so that unwanted vibrations resulting therefrom are substantially 
suppressed. This technique has at least two advantages. First, it is relatively 
easy to implement from a computational point of view because the problem is 
separated into two problems, namely generating a good input for a simple system 
and separately generating a good Shaper™ to handle flexible, and possibly some 
rigid, modes. Second, solutions obtained by separately Shaping''''^ the input are 
more insensitive to system parameter variations. 

Shaping of the input may be performed in a maimer described 
above using, e.g., an Input Shaper™. a filter, etc. Accordingly, a detailed 
description of shaping the input is omitted here for the sake of brevity. With 
regard to generating a good input, since voltage levels for the input are already 
known (i.e., the voltage levels are either full positive and full negative), the only 
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variables that need to be determined are switch times for the voltage input. 
Switch times can be determined using a partial fraction expansion model of a 
system, such as equations 19, in which F,(t), FjCt) and F^ii) correspond to step 
responses to a unit voltage in each mode of the system. For a third order model 
such as this, three transitions between full positive and full negative are required. 
Switch times, namely ty, t, and t^nd, for the transitions are determined as follows. 
Using equations 14, reproduced below as 39, 

^^2i'eJ = m^r,, - to) - 2KoF,(r„, - f,) - 2V^^(t^„, - g (39) 



and constraint equations 15, reproduced below as 40, 

XF,{U - 0 (40) 



it is possible to solve for the three unknown transition times t^, and t^^^. A 
change of variables, i.e., 







-^0 




et^ 


~ *end 




(41) 


et. 









fiirtlier simplifies the solution for the switch times. In this case, eti corresponds 
to an elapsed time from onset of a first step until the end of the move, et, 
corresponds to the elapsed time from onset of the second step until the end of the 
move, and etj corresponds to the elapsed tune from onset of the third step until 
the end of the move. 

In view of die foregoing, equations 39 reduce to the following: 
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0 = V,F,(et,) - 2V,F,iet,) ^ 2V,F,{et,) (42) 
0 = Wet,) - 2V,F,(et,) ^ 2V^F,iet,) . 

Since equations for Fi(t), F2(t) and F3(t) contain exponential terms, a numerical 
solution to the above equations for XFi(t^d)- ^^li^^J and XF3(t,J is not readily 
available. However, certain simplifications can be made to the equations in order 
to provide an iterative numerical solution, and thus to allow them to be solved for 
5 transition times tj, t, and tg^d in real-time. More specifically, the duration of the 
third step of the input to a third order system does not vary significantly for 
relatively long moves. Therefore, for' long moves, the third step duration is 
nearly independent of values of etj and et,. Accordingly, the second equation (for 
XF2(te„u)) can be solved for etj using guesses for etj and et,. The value of etj 

10 obtained in this manner is a close approximation to its exact value. The first 
equation (for XFi(tend)) is used to eliminate et, from the third equation (for 
XFsCtgJ), whereafter the third equation is solved for et jas a function of ei .3 The 
value obtained for et, is also a close approximation of its exact value. The first 
equation is then used to calculate an approximation for etj based on the 

15 approxinaate values of et, and etj. If more accuracy is desired, another pass may 
be made through the equations. That is, the approximations calculated for et, and 
etj can be substituted back into the foregoing equations, and the process repeated 
for several iterations, where each iteration provides an increase in accuracy. In 
practice, four iterations through the equations resulted in 16 digit accuracy for et,, 

20 et2 and ti^. Once these values are determmed, values for the transition times tj, 
t2 and tend can be calculated using the equations shown above. 

The iterative solution presented above for determining switch times 
will converge faster for better initial guesses of et, and et,. On short moves 
especially, convergence can be slow if initial guesses for et, and et2 are not good. 

25 Accordingly » there exists a need for a way to provide a good initial guess. One 
such technique is to tabulate solutions for a range of short moves. Such a table 
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may be then used to obtain a good guess, e.g., by looking-up a guess or by 
interpolation. Alternatively, a curve can be fit to solutions for a range of short 
moves, and initial guesses determined based on the curve. For example, it has 
been determined experimentally that a first switch time is a monotonically 
increasing function of move distance. Curves for other switch times can be 
determined in a similar manner. 

The foregoing techniques for determining switch times can be 
expanded to include parameters other tlian just movement distance. For example, 
if resistance change is to be included in the solution, an initial guess for a switch 
time value can be determined based on both movement distance and resistance. 
Any number of system parameters can be used in this manner to determine an 
initial guess. In such cases, the initial guess can be determined based on a 
multidimensional table which includes these parameters, or computed using an 
empirically generated function. Switch times can also be generated by taking into 
account system parameters in addition to those set forth above. In such cases, 
conservative parameters may be used so that the actual system will move faster 
flian the theoretical model on which the switch times are based. That is, in this 
case, the feedback will scale back the output voltage and the system will almost 
never be commanded past saturation. 

As an alternative to the above approaches which rely on an initial 
guess, short move solutions can be determined based on the partial fraction 
expansion equations. Specifically, the exponential terms of the partial fraction 
expansion equations can be expressed in series form. The equations can then be 
solved for the short move solutions to the switch time problem. 

A hardware implementation of this embodiment of the invention 
may be dictated by the design of many disk drive controllers. Specifically, in 
such disk drive controllers, the update rate of a processor in the controller is not 
fast enough to output voltages at precise desired time intervals. Accordingly, a 
hardware solution may be implemented, in which each next switch time is loaded 
into a hardware counter chip and a high frequency clock counts intervals until the 
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desired switch time is achieved. At such time, the output of the counter chip 
changes and the new voltage value is output to the disk drive system. 

Still another alternative to the forgoing method for determining 
switch times for a system input is provided in U.S. Provisional Application No, 
60/076,951, entitled "Creating Time-Optimal Commands For Linear Systems", 
filed on March 5, 1998 and its corresponding national application filed 
concurrently herewith. The contents of these two applications are hereby 
incorporated by reference into the subject application as if set forth in full. 

7.0 Shaping"'^'^^ Near or Bevond Nvquist 

As noted above, control techniques in which trajectory is specified 
as a function of time involve determining an input which will generate a good 
move. As also noted above, this input may be determined using any other 
techniques, such as Input Shaping"^ filtering, etc. In addition to these 
techniques, the present invention also provides a technique for shaping inputs 
having frequencies near or beyond the nyquist frequency. Shapers™ used to 
perform this shaping are referred to herein as "NBN" Shapers'^'^. 

Systems that use discrete time or digital control face a fundamental 
limitation due to the effects of sampling. That is, the sampling theorem states 
that any frequencies greater than the sampling frequency (or Nyquist frequency) 
will be aliased to other frequencies. Therefore system feedback control carmot 
suppress vibrations that are at a frequency greater than the Nyquist frequency. 
Because Input Shaping^^* is a feedforward technique, it does not have this same 
limitation. Input Shapers"™ can be designed to suppress vibrational frequencies 
that are near, or even well above, the Nyquist frequency. This section will 
demonstrate how to design Input Shapers*'^'"^* which limit vibrational frequencies 
that are near or above the Nyquist frequency. 

As background, disk drives systems with embedded servo 
information have a maximum position update rate that is determined by the 
density of the servo information and the drive's rotation speed. To demonstrate 
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the difficulty in solving for an Input Shaper"^^, an example drive is presented. 
Assuming that the servo information on such a disk drive is updated at 5 KHz, 
and the drive's first resonant mode is at 2 KHz, only 2-f- servo cycles exist for 
each vibrational cycle. This presents a computational problem. That is, only 2 
5 unknowns exist per servo cycle and there are many more constraints to be solved 
(at a minimum no vibration constraints, among others). Accordingly, the 
problem is over-constrained. 

However » more than one cycle can be used to solve for a system 
input. By expanding the number of unknowns until the problem becomes 

10 ' solvable, a solution can be obtained. Solutions obtained are unusual because they 
take more than one cycle of vibration to reduce vibrations to an acceptable level. 
Additionally, careful attention must be paid to insensitivity constraints because 
NBN solutions will be less sensitive than conventional solutions since fewer 
cycles of vibration are allowed to occur before cancellation. An alternative to 

15 this approach is to use an output or feedforward command rate that is different 
from the servo rate. This technique is presented in section 7.1 below. 

The basic principle of the NBN shaper can be understood by 
examining the design of a three pulse shaper for a continuous system with zero 
damping. If we assume that the system is continuous, the shaper unpulses can 

20 occur at exact times. When the damping ratio is zero, the araplimdes of the 
impulses for a three pulse shaper are [.25 .50 .25]; the first pulse occurs at time 
zero, the second pulse at one half period of vibration, and the third pulse at one 
period of vibration. Figure 32 shows a plot of the impulses of the three pulse 
shaper along with a sine wave representative of the undamped vibration. The 

25 combined plot shows ±e phasing of the impulses relative to the sine wave. In the 
continuous case, the pulses occur at the times shown hi Figure 32. In the digital 
case, the pulses are constrained to occur at integer multiples of the sample period, 
i.e. at tlie sample tunes. If the sample period is much shorter than the vibration 
period, then each of the shaper pulses can be split between the nearest sample 

30 limes to achieve the same effect. 
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The Nyquist rate is equal to one half the sample rate and represents 
the lowest frequency that can be sampled by a digital system without alias. Thus, 
vibration occurring at the Nyquist rate is just at the theoretical limit of what input 
can be processed by a digital system without alias, and results in two sample 
5 points per sine wave period. Two sample points per period is exactly die ideal 
spacing for a three pulse shaper. Figure 33 shows a plot of the digital pulse 
spacmg for a three pulse shaper where the digital sample rate (4 KHz) results in 
the Nyquist rate (sample_rate/2 = 2 KHz) being equal to the vibration frequency 
{2 KHz). Thus, the three pulse shaper is easily unplemented to attenuate 

10 frequencies at the Nyquist rate. 

It is well known that vibration with frequency in excess of the 
Nyquist rate caimot be sampled by a digital system without aliasing, and that 
therefore a closed-loop digital sampled data system cannot use closed-loop control 
techniques to attenuate vibration at high (beyond Nyquist rate) frequency. 

15 However, the open-loop (feedforward) shaping technique can be applied to 
attenuate vibration at frequency well beyond the Nyquist rate. Shown in Figure 
34 is the pulse spacing for a three pulse shaper and a sine wave representative of 
vibration. The combined plot shows the timing of the shaper pulses relative to 
the phasing of the sine wave. The frequency of vibration is again 2 KHz, the 

20 sample rate for this example is 1.333 KHz, and the associated Nyquist rate is 
1.333KHz/2 = 0.666KHZ, The frequency of vibration is well in excess of the 
Nyquist rate (2KHz/0.666KHz = 3) for this case. 

Examination of the phasing of the sine wave relative to the timing 
of the shaper pulses reveals the same relation for the sub-Nyquist rate example 

25 as for the Nyquist rate example. That is, the first shaper pulse occurs at 0**, the 
second shaper pulse occurs at 180** and at each 360"* thereafter, and the third 
shaper pulse occurs at 360° and at each 360*' thereafter. Since the sine wave 
phase repeats each 360°, the superposition of responses yields the same result in 
both cases. It should be noted that the sub-Nyquist rate shaper will be more 

30 sensuive to errors in the vibration frequency estimate. The increased sensitivity 
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is due to the fact that as more periods of vibration elapse between pulses, larger 
errors in phase (for a given error in frequency) will accumulate. 

7. 1 Effect of the Rate of Command Output 
5 The rate that is chosen for issuing digital commands to the system 

is extremely important. A rate should be chosen that makes sense for the 
panicular system being considered. In some cases, this rate is simply the servo 
rate of the control system. However, for disk drives, this rate must usually be 
higher because of the restrictions imposed by the use of embedded servo 

10 information for control. For systems in which the highest frequency of vibration 
to be suppressed is considerably lower than the servo rate, the servo rate is 
adequate for use as an output rate. For systems such as those considered in 
section 7.0 (such as most disk drives), the rate must be higher. 

One approach is to send output commands to the system at a 

15 multiple of the servo rate. In this formulation, the output to the physical system 
might be performed at two, three, or more times the rate of the servo. In our 
example system, the servo rate is SKHz. The highest mode to be suppressed is 
2KHz. If the output to the system is performed at triple the servo rate (ISICHz), 
there are over seven digital outputs per vibration cycle; therefore, the shaping 

20 formulations become simple. Several possible architectures can be considered. 
Two possibilities are presented below. 

In the first example, there is a separate output path from the 
controlling computer chip to the hardware for feedforward purposes. An example 
of this is a disk drive in which a voltage command can be fed forward directly 

25 to the amplifier while current commands are output by the closed-loop servo. In 
order to actually implement this, the following steps must be performed. At the 
servo rate (5KHz). these steps comprise reading the sensors, executing the servo 
calculations, computing the servo output, and outputting this value to the servo 
output channel. At the higher rate (15KHz). these steps comprise computing or 

30 looking-up the next trajectory value, executing the shaping method (if the 
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correction for unwanted dynamics is not already included in the trajectory), and 
outpuuing the feedforward term to the feedforward output channel. 

In the second example, there is only a single output path from the 
controlling computer chip to the hardware. An example of this is a disk drive in 

5 which a current is the output generated by the closed-loop servo and current is 
also used in the feedforward term. In order to actually implement this, the 
following steps are performed. At the servo rate (5KHz), these steps comprise 
reading the sensors, executing the servo calculations, computing the servo output, 
and saving this value for use by the section below. At the higher rate (15KHz), 

10 these steps comprise computing or looking up the next trajectory value, executing 
the shaping method (if the correction for unwanted dynamics is not already 
included in the trajectory), and adding the feedforward term to the feedforward 
term to the servo output computed above. 

There are many related configurations which either combine the 

15 above two approaches or are variations of these approaches. One such variation 
uses continuous or near continuous time outputs from the feedforward term. If 
the control hardware for a disk drive or any other machine is constructed so that 
outputs can be scheduled in time, the feedforward outputs can run essentially 
asynchronously from the servo loop. In one such implementation, the voltages 

20 to the disk drive might be updated by a circuit at a high clock rate (e.g., 2 MHz). 

By loading a counter register with the time of the next voltage output, the 
hardware circuit can count down until the appropriate time and output the next 
voltage value. This enables an optimized trajectory to be computed using any of 
the techniques herein (or by any other process), and to be output at a rate 

25 sufficiently high so that the system can be commanded to move in an optimized 
(or near optimized) fashion without effecting the servo rate. 

An additional variation on the approaches presented herein is to use 
one of the shaped switched formulations so that only full scale voltages (or any 
other parameter) are output. For example, it is possible to use unity magnimde 

30 Shapers*^*"^* to take a switching voltage command signal designed only for a rigid 
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body system, and to modify the switching vokage command signal so that the 
switching voltage command signal moves the system with little vibration. This 
approach enables the circuit to become a simple counter driven switch that 
changes its output at a preset time. 

A further variation is to use an interrupt-driven approach in which 
times to output values are scheduled and an interrupt is generated (using state of 
the art computer techniques) at the appropriate time. At the time of the interrupt, 
an appropriate output (or voltage switch) is provided by the controller chip or 
other device. Single or multiple hardware timers may be used to interrupt a DSP 
or microcontroller in the disk drive. Still another variation is to use a separate 
device, such as a programmable logic array, that is issued a single command 
which, in turn, generates a series of precisely -timed commands to the current 
loop. The DSP or microcontroller can be progranmied in an asynchronous 
control loop that maintains a list of switch times and required transitions, 

8.0 PV Approach 

A PV table comprises a trajectory curve which defines the velocity 
of a system component, such as head 76, in terms of that component's position, 
i.e., the distance that the component has yet to go in order to reach its final 
position. To use a PV table, an actual position of the component is compared to 
its desired position. The difference between these two values is then used to look 
up a velocity in the PV table for the component. Many disk drive systems rely 
on a PV table to determine their inputs. These inputs, however, can result in 
unwanted vibrations in the system. The techniques presented below can be used 
to reduce and/or limit such vibrations. 

This aspect of the invention is a method of generating a trajectory 
for inclusion in a PV table which is used m controlling movement of a component 
in a dynamic system. The method generates a trajectory in accordance with a 
technique for reducing unwanted dynamics of the component, i.e., generating a 
good move. The trajectory is then stored in the PV table, and the component is 



-63- 



wo 99/45535 



PCT/US99/04910 



controlled to move along the trajectory stored in the PV table. 

In more detail, a good move is generated for the system using any 
of the techniques described in section 5 above. Specifically, the system is 
modeled, a "good" input is applied to the system, and the system's reaction to the 
5 good input is stored in a PV table as a trajectory for the system. In this mamier, 
good moves can be re-cast into a PV table for a system, thereby providing a way 
to reduce unwanted system dynamics while staying within a PV formulation. 
Taking this approach one step further, a series of good moves can be generated 
using this technique, and the series of good moves can be stored in a series of 
10 parameterized PV tables. 

9.0 PV Feedback Approach 

In prior art disk drive systems, the velocity determined from a PV 
table is merely fed forward into the system, as shown, for example, in the 

15 conventional disk drive system diagram depicted in Figure 13. This could lead 
to unwanted system vibrations. This aspect of the invention, however, shapes the 
velocity before it is provided to the system and, in addition, performs inverse 
shaping m the system's feedback loop, thereby correcting the feedback so as not 
to co-opt the operation of the PV table. 

20 Figure 14 shows an example of a system diagram for a disk drive 

that is controlled in accordance with this aspect of the invention. The disk drive 
includes PV table 84, shaper 86, and inverse shaper 87. Shaper 86 is placed in 
the feedforward loop following PV table 84 so as to shape the velocity output 
therefrom. Inverse shaper 87, which comprises the inverse of shaper 86, is 

25 placed in the feedback loop. The relationship between inverse shape 87 and 
shaper 86, in either the time domain or the frequency domain, is as follows: 

if shaper - A/B, then shaper^ = B/A. 

30 The combination of shaper 86 and inverse shaper 87 reduce unwanted system 
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dynamics, such as vibrations of head 76. Specifically, shaper 86 alters the 
velocity input to the system so as to reduce vibrations in head 76 resulting from 
system movement, and inverse shaper 87 operates to restore loop stability to what 
it was prior to insertion of the shaper 86 in the forward loop. Shaper 86 and 
inverse shaper 87 are not limited to the Input Shapers*^*^ and inverse Input 
Shapers"^^^ described in U.S. Patents Nos. 4,916,635 and 5,638,635, In fact, 
shaper 86 and inverse shaper 87 can be any pole-zero dynamic block, IIR filter, 
FIR filter, tapped continuous delay line filter, or the like. 

Figures 15 and 16 show relationships between inputs and outputs 
of a shaper ("IS'*) and inverse shaper ("IS*'"), respectively, formed using an FIR 
filter, with A as the input. Specifically, relationships between the inputs and 
outputs of IS and IS'' are as follows: 



yn 



N 
/=0 



(43) 



Figures 17 and 18, show relationships between inputs and outputs 
of a shaper and an inverse shaper ("IS"'"), respectively, which are formed using 
a tapped continuous delay line filter. Specifically, relationships between inputs 
and outputs of IS and IS ' are as follows: 



(44) 



Figures 19 and 20, show relationships between inputs and outputs 
of a shaper ("IS") and an inverse shaper ("IS'*'*), respectively, which are formed 
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10 



15 



using an IIR filter. Specifically, relationships between inputs and outputs of IS 
and IS*^ are as follows: 



and corresponding apparatus for using a PV table to control a component, such 
as head 76, in a dynamic system, such as disk drive 10. The invention generates 
a position variable for the component and determines a velocity command for the 
component using the PV table, where the velocity conmiand is determined based 
on the position variable. The velocity conmiand is then shaped using an Input 
Shaper^'*^ or the like in order to generate a shaped velocity command, whereafter 
the component is controlled to move based on the shaped velocity conunand. In 
the invention, a preset position of the component is compared to a measured 
position of the component in order to determine the position variable, and inverse 
shaping is performed on the measured position prior to comparing the measured 
position to the preset position. 

10.0 N-Dimensional PV Table 

The present invention also provides an enhanced PV table which 
may, or may not, be used in connection with the foregoing feedback approach and 
the foregoing techniques for generating good moves. As background, the same 
PV table may not produce accurate results for different distance moves, among 
other things. To address this problem, conventional disk drive systems typically 
include one, two or maybe even three PV tables, i.e., one table for short moves, 
one table for medium-distance moves, and one table for long moves. Which table 
is selected depends upon the distance of the desired move. Since only three 



1 ^ 




(45) 




To summarize, this aspect of the invention is directed to a method 
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discrete values are provided for the movement distance, the system may not reach 
its final destination as quickly as desired. 

The present invention provides a belter approach, namely to 
generate an N-dimensional (N>2) PV table in which velocity is a function of 
position and a function of a third (fourth, fifth, etc.) variable, such as movement 
distance, resistance, temperature, etc. This approach effectively provides 
increased movement speed, since the requested velocity values more closely 
represent physical limitations of the disk drive for substantially all movement 
distances. Moreover, this approach makes more sense than the conventional 
formulation, especially in the disk drive area, since disk drive systems are third 
order in nature and a conventional PV approach is second order in nature. 
Although the present invention may use a PV table having any number N (N > 2) 
of dimensions, for the sake of simplicity, the examples provided herein refer 
usually to a three-dimensional PV table. 

In order to determine trajectories for a PV table, a dynamic model 
of the system must be determined. An example of such a model is partial 
fraction expansion equations 19. A system input, which is a function of time and 
which will move the model a desired distance, is then generated based on the 
model. For example, the system input may consist of a series of switch times 
and voltage levels, or a digital sequence of voltage levels. The input, in 
conjunction with the system model, is used to simulate the desired motion of the 
system. Any form of the model that yields both a position and a velocity output 
is suitable. As an example, the third order voltage-to-position transfer function 

His) = ^ * * (46) 

could be transformed and manipulated so that position and velocity are available 
as outputs. Using the commutative property of multiplication, the transfer 
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function can be rearranged and then broken up to yield a "block diagram" with 
voltage input and position and velocity outputs as follows: 

Voltage - ^ — - — - — - Velocity 

J + ^7 s + b ^4*7^ 

Velocity -* — - Position 
s 

The above "block diagram" is then used as the basis of a computer 
simularion that yields the velocity and position outputs corresponding the voltage 
input. The integration is performed forward in time. If both the form of the 
input (e.g. three switch times and three voltage levels) and the model are simple, 
then direct analytical solution of the system model to time domain equations can 
be used as an alternative approach to simulation. The velocity, and position 
outputs of the simulation are then tabulated to form one version of the PV table. 

The present invention enables three approaches to PV table 
generation. First, a shaped input can be integrated forward to generate a PV 
table - in the above block diagram a shaper would be inserted following the 
voltage input. Second, an unshaped input can be integrated forward and an 
output of the integration can be shaped and stored in a table. Third, a vibration- 
reducing trajectory can be computed using a constraint-based approach, this input 
can be integrated forward, and the output thereof stored in a table. Prior to the 
derivation of the partial fraction expansion equations shown above, there was no 
way to generate an input that was guaranteed to get the system to its final 
destination; therefore forward integration was not a viable approach. Reverse 
integration from a final condition was commonly used. 

Using the generated input, system outputs are then determined over 
lime so as to generate values of position as a function of time and velocity as a 
function of time. These values are then re-cast as velocity as a function of 
position. To generate a three-dimensional PV table, the above process is repeated 
for the third variable in order to generate velocity as a function of both position 
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and the third variable. The process may then be repeated for additional variables. 

Since a three-dimensional PV table may require large amounts of 
memory space, and therefore may not be practical in all circumstances, the 
present invention provides several ways in which to approximate operation of a 
5 three-dimensional PV table using a series of PV trajectories. For example, the 
invention generates a series of PV trajectories for various movement distances, 
stores these trajectories, and selects trajectories corresponding to a desired 
movement distance while ignoring all other trajectories. Thus, in this case, one 
of the series trajectories is selected from the position-velocity table based on the 

10 desired movement distance and the system is controlled accordingly. 

Another way to approximate a three-dimensional PV table is to 
perform curve fining on a series uajectories in a PV table, such as those shown 
in the PV table of Figure 21. With regard to Figure 21, it is noted that, in the 
reahn of computer disk drives, an actual target velocity is typically not included 

15 in the PV table. Instead, a velocity command, comprising a target velocity plus 
an additional tenn relating to the actuator arm's effort to stay on track, is 
included in the PV table. With regard to curve fitting, since the curves shown 
in Figure 21 are substantially regular, it is possible to parameterize the curves and 
to generate a function which can be used to generate a single PV trajectory to 

20 control the system. 

Figure 22 shows examples of parameterized curves which were 
generated from the curves shown in Figure 21. For the example shown in Figure 
22, an estimate of curves between the lowest (3000) and the highest (100) curve 
could be obtained by taking a weighted average of the lowest and highest curves 

25 based on movement distance. Consequently, in this example, an entire series of 
curves can be estimated by taking a weighted average of only two curves and 
using original movement distance to adjust the weights in the weighted average. 
Another way to approximate a three-dimensional PV table is to generate a custom 
PV table for each movement distance in real-time. This can be done using the 

30 system equations 19 for a computer disk drive set forth in section 4 above. 
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It is also possible to modify a series of trajectories in a PV table 
based on measurements of system parameter variations, such as variations in 
temperature, resistance, etc. This is done by estimating system parameters while 
motion is in progress, and determining v/hether the system parameters have varied 

5 from predetermined system parameters. Trajectories in the PV table are then 
modified based on any determined system parameter variations, whereafter the 
modified trajectories are stored in the PV table. Any method of estimating the 
system parameters may be used in connection with the present invention, 
including a conventional parameter estimator. For example, as a disk drive heats 

10 up, resistance therein changes. By looking at prior move responses, it is possible 
to estimate the resistance, to determine whether the resistance has changed, and 
to recalculate the drive's PV table, if necessary. 

More specifically, m a disk drive system, during a first part of a 
move, the system is in saturation and, therefore, the PV table is inactive. This 

15 corresponds to issuing a step command in voltage to accelerate movement as 
rapidly as possible. During this acceleration phase, the acmal current flowing 
through a voice coil motor used to convert electrical energy to mechanical energy 
can be measured. Knowing this, it is possible to perform an energy balance 
calculation to determine energy losses in the system due to resistance. The 

20 energy lost due to resistance equals the energy input to the disk drive minus 
kinetic energy of the voice coil motor minus energy stored in the voice coil 
mductor. This determination is made assuming that the inertia of the voice coil 
motor and the inductance of the voice coil are known or can be determined. 

In real-time, the following steps are performed in order to estimate 

25 resistance: (i) for each time step, measure the acmal current and the actual 
voltage applied to the drive's voice coil motor, (ii) integrate power in the drive 
(determined by multiplying voltage by current) over time to determine the total 
energy input to the system, (iii) calculate the instantaneous kinetic energy of the 
voice coil motor, (iv) calculate the instantaneous energy in the motor inductance, 

30 (v) calculate the integral of drive current squared over time, (vi) determine the 
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estimated resistance by subtracting, from the energy input, the kinetic energy and 
the energy in the inductor divided by the integral of the current over time. Thus, 
after each integration, an estimated resistance for the disk drive is calculated. 
This can then be used to update the PV table. 
5 Similar energy of momentum balance equations can be used to 

estimate other system parameters. For example, the following equation can be 
used to estimate motor torque: 



K, = — . (48) 



idt 



This equation can be solved in real-time, and can be solved over a whole move 

10 or incrementally starting at any lime during the move. 

The foregoing techniques for generating PV tables can be 
calculated off-line, thereby making it possible to parameterize or curve-fit PV 
tables as a function of system parameters as well as movement distance. 
Consequently, it is then possible to retrieve a PV table based on movement 

15 distance and system parameters. Moreover, it is noted that an N-dimensional PV 
table can be used outside of the context of the vibration reduction methods 
described herein, and that a non-dimensional PV table can be used in connection 
with all aspects of the present invention. A non-dimensional PV table comprises 
a PV table in which the independent variable on the X-axis of the PV table has 

20 been scaled using a predetermined scaling factor. Such a factor may be, for 
example, a maximum number of tracks moved, or the like. 

11.0 Shape on Transient ("SHOT") 

Conventional implementations of the PV approach to controlling 
25 a disk drive encompass variations of the following steps: (i) ramping the current 
command to saturation (various techniques are used for this); (ii) commanding 
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full, saturated ctnrent until some braking condition is satisfied, (iii) transitioning 
to full, saturated deceleration until the PV curve is crossed; (iv) following tlie PV 
curve until the drive is close to the final destination; and (v) turning on a terminal 
controller to bring the system on track. The preponderance of vibration, both 

5 mechanical and acoustic, is created by the sudden transitions at the beginning and 
middle of the PV move. Consequently, Input Shaping™ (or any other vibration- 
reducing technique) need be applied to only these portions of the trajectory. 

A new sequence of operations, using SHOT, is as follows: (i) 
Input Shape™ the current command to saturation; (ii) command full, samrated 

10 current until some braking condition is satisfied; (iii) Input Shape™ the transition 
to full, saturated deceleration until the PV curve is crossed; (iv) follow the PV 
curve until the drive is close to its final destination; and (v) turn on a terminal 
controller to bring the system on track. Since the transitions to saturation are 
both just steps, the shaped transitions become shaped steps (or staircases). 

15 Shaping can be rendered less effective if the input results in 

saturation and the command generator is unaware of the saturation. There are 
several methods for dealing with this situation. It is possible to command the 
system into saturation using estimates of when saturation should occur based on 
previous moves. For example, by estimating the maximum voltage and the coil 

20 resistance (or the quotient of maximum voltage/coil resistance) the maximum 
current in saturation can be obtained. Once the system is shaped into saturation, 
the saturation current can be calculated by measuring the acceleration and the 
velocity, since current is related to acceleration by a constant and the effective 
voltage is die saturation voltage minus the back EMF constant times the velocity. 

25 Based on these calculations, it is possible to compute the actual 

saturation level. It is then important to command only enough of an input to keep 
the system in saturation and not exceed this level by too great of a value. If this 
value is exceeded by a considerable amount, the shaping that is applied when the 
command is reversed to slow the system down will be "chopped-off" at the 

30 beginning. For example, if an output command of 5 volts from a DAC will 
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saturate the current amplifier in steady state at the current velocity, then a 
command of 6 volts would be undesirable. During the transition to deceleration, 
if this error is left uncorrected, the shaper will reduce the voltage to, e.g., 5.5 
volts as the first part of the shaping process, yet the amplifier will remain 

5 saturated. Therefore, the effect of shaping is degraded. 

While the SHOT technique is well suited for adding to the PV 
approach, it is by no means limited to PV. Any method of commanding a disk 
drive to move that causes significant transitions in the command to the system is 
a candidate for SHOT. The user need simply use shaped transitions instead of 

10 sudden or ramped transitions. The key to this approach is that while sudden 
transitions cause vibrations and noise, and smooth transitions often cause 
vibrations and noise (although less so than sudden transitions); shaped transitions 
can be constructed to yield any level of vibrations and noise. This makes the use 
of shaped uansitions desirable because the effects are predictable and the move 

15 is lengthened by the minimum amount. SHOT can also be used for transitions 
caused by sudden gain changes as well as command changes. In this application, 
the gain is altered in a shaped manner from one level to another using any one 
of the shaping techniques described herein or from the literature. 

20 12.0 Conclusion 

At this point, it is reiterated that all aspects of the present invention 
can be implemented in hardware, software or a combination of thereof. 
Moreover, although the invention has been described in the context of computer 
disk drives, it may apply equally to other types of data storage devices (of which 

25 a computer disk drive is one), including, but not limited to, optical drives, tape 
drives, dual-actuated disk drives, and holographic storage devices which read 
from, and write to, data storage media other than magnetic disks. 
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The present invention has been described with respect to particular 
illustrative embodiments. It is to be understood that the invention is not limited 
to the above-described embodiments and modifications thereto, and that various 
changes and modifications may be made by those of ordinary skill in the art 
without departing from the spirit and scope of the appended clauns. 



-74- 



wo 99/45535 PCT/US99/0491 0 



APPENDIX A 

% BEGINNING OF PSEUDO CODE 

5 % compute scale factor A, and time constants a, b from physical system 
% parameters 

A = Vmax * Kt / (Re * Rm + Kt * Kb) * 1 * k; 

10 pi = 1/Jm/Ie * (-le * Rm - Re * Jm + sqn(Ie"2 * Rm"2 - 2 * Re * Rm * le 

* Jm + Re"2 * Jm"2 - 4 * Kt * Kb * le * Jm)) / 2; 

p2 = 1/Jm/Ie * (-le * Rm - Re * Jm - sqrt(Ie"2 * Rm"2 - 2 * Re * Rm * le 

* Jm + Re"2 * Jm"2 - 4 * Kt * Kb * le * Jm)) / 2; 

15 a = max(-pl,-p2) 
b = min(-pl,'-p2) 

% make initial guesses for step durations 

20 etl = I; 

et2 = .005; 
et3 = I; 

% set maximum iteration count 

25 

Nmax = 1000; 

for j = l:Nmax 

Sfc save old values of step lime imervals 
30 etSold = et3; 
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et2old = et2; 
etlold = etl; 

% iterate for switch times using fixed voltage level Vmax 

et3 = -log(l,0 / 2.0 • exp(-etl * a) / 2 + exp(-et2 * a)) / a; 

et2 = 1/b * log(2.0) + 3 * et3 - 1/b * Iog(2 * exp(l/A * b * X) * 

exp(et3 * b) - sqrt( 4.0) * sqrt(exp(l/A * b * X)) * exp(et3 * b) 
* sqrt(exp(l/A * b * X)+exp(et3 * b)"2 - 2 * exp(et3 * b))); 

etl = - (-2 * A * et2 + 2 * A * et3 - X) / A; 

if norm((et3old - etS et2old - et2 etlold - etl], inf) < = eps * 2 
break 

end 

if j==Nmax 

error([ 'error - failure to converge after num2scr(Nmax)/ 
iterations']) 

end 

end 

% round up pulse duration to nearest sample interval, 

% convert to intervals between steps to make sure that voltage 

% requirements will not increase (beyond Vmax). 

dtl=ceil((etl - et2) / dt) * dt; 
dt2=ceil((et2 - et3) / dt) * dt; 
dt3=ceil((et3) / dt) * dt; 

etl23 = [etl, et2, et3] 

% convert back to total step duration. 
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ell = dtl + di2 + dt3; 
ei2 = dt2 + di3; 
ei3 = dt3; 

% In the following, the original constraints equations involving XFl. XF2, 
% and XF3 have been modified to include a variable voltage level applied at 
% each step (instead of the fixed maximum (+/-) Vmax). 

% The original equations for XFl, XF2. and XF3 follow: 
% XF,(t„,) = VoF,(t,„, - g - 2V,FXt,„, - 1.) + 2VoF,(t^, - t^) 
% XF^d,^) - VoF,(t,,,, - g - 2VoF2(t,„, - t,) + 2VoF,(ten, - t.) 
% XF3(t,^) = VoF3(t,^ ~ to) - 2VoF,(te„a - t^) + 2VJP,{i,,, - 

% And the modified equation including adjustable relative levels of voltage 
% LI, L2 and L3 are: 

% XF,(t,J = L»VoF,(t,„, - tc) - L,VcF It^,- t ) + L 3V f (i t ) 
% XF,(teJ - LiVoF,(t,^ - to) - UVoF,(t,„, - t,) + L3VoF,(t,,, - t,) 
% XFjd^) = L,VoF3(t,,^ - to) - UVoF^Ct^, - t,) H- LsVoFiCt,^ - t,) 

% And the corresponding constraint equations are: 
% XF,(t,^) = Finalpos 
% XF,(te„,) = 0 
?c XF,(i,J = 0 

% Where all of the times indicated have discrete values, e.g. corresponding to 
% the controller update rate. 

% It should be noted that after the digital switch times are fixed, the constraint 
% equations derived from the equations above form a linear set of equations in 
% the unknown relative voltage levels LI, L2 and L3 and any standard linear 
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% meihod can be used to solve for the relative voltage levels, in the equations 
% for (LI, L2 and L3) that follow, the solution was obtained by algebraic 
% means (and are not panicularly compact.) 

% compute new relative voltage step levels 

% LI, L2 and L3 are nominally assigned to "1". "-2" and "+2", respectively 

si = X * (exp(-et3 * b) * exp(.et2 * a) + exp(-et3 * a) -h exp(-et2 *b ) - 
exp(-et2 *b) * exp(-et3 * a) - exp(-et2 * a) - exp(-et3 *b )); 

s2 = 1 / (et2 * exp(-etl * b) * exp(-et3 * a) -f exp(-et2 * b) * et3 * 
exp(-etl * a)- ct2 * exp(-et3 * a) - et2 * exp(-etl * b) -et3 * 
exp(-etl * a) - exp(-et2 * b) * et3+exp(-et3 * b) * etl * exp(-et2 * a) 
+ exp(-et3 * a) * etl + exp(-et2 * b) * etl - exp(-et2 * b) * etl * 
exp(-et3 * a) - et3 * exp(-etl * b) * exp(-et2 * a) - exp(-et2 * a) * etl - 
exp(-et3 * b) * etl - exp(-et3 * b) * et2 * exp(-etl * a) -f et3 * 
exp(-etl * b) 4- et2 * exp(-etl * a) + exp(>et3 * b) * et2 + et3 * 
exp(-et2 * a)) / A; 

LI = si * s2; - 

si = 1 / (et2 * exp(-etl * b) * exp(-et3 * a) + exp(-et2 * b) * et3 * 
exp(-etl * a) - et2 * exp(-et3 * a) - ei2 * exp(-etl * b) - et3 * 
exp(>etl * a) - exp(-et2 * b) * et3 -f exp(-^et3 * b) * etl * 
exp(-et2 * a) + exp(-et3 * a) * etl + exp(-et2 * b) * etl - 
exp(-et2 * b) * etl * exp(-et3 * a) - et3 exp(-etl * b) * 
exp(-et2 * a) - exp(>et2 * a) * eil - exp(-et3 * b) * etl - exp(-et3 * b) * 
et2 exp(-etl * a) + et3 * exp("ecl * b) + et2 * exp(-eil a) + 
exp(-et3 *b ) * et2 + et3 * exp(-et2 * a)) * X; 

s2 = (exp(-et2 * b) * exp(-eil * a) - exp(-etl * a) - exp(-et2 * b) - 

exp('eil b) * exp(-et2 * a) -f exp(-etl * b) + exp(-et2 * a)) / A; 
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L3 = sl*s2; 

si = exp(-etl * a) - exp(-et3 * a) + exp(-et3 * b) - exp(-etl * b) - 
exp(-et3 * b) * exp(-etl * a) + exp(-etl * b) * exp(-et3 * a); 

s2 = X / (et2 * exp(-etl * b) * exp(-et3 * a) + exp(-et2 * b) * et3 * 
exp(-etl * a) - et2 * exp(-et3 * a) - et2 * exp(-etl * b) - et3 * 
exp(-etl * a) - exp(-et2 * b) * etS + exp(-et3 * b) * etl * exp(-ec2 * a) 
+ exp(-et3 * a) * etl + exp(-et2 * b) * etl - exp(-et2 * b) * etl * 
exp(-et3 * a) - et3 * exp(-etl * b) * exp(-et2 *a ) - exp(-et2 * a) * etl- 
exp(-et3 * b) * etl - exp(-et3 * b) * et2 * exp(-etl * a) + et3 * 
€xp(-etl * b) + et2 * exp{-etl * a) + exp(-et3 * b) * et2 + et3 * 
exp(-et2 * a)) / A; 

L2 = si * s2; 

% convert accumulated voltage steps to sequential voltage level 

VI = Vmax * (LI); 

V2 = Vmax * (LI + L2); 

V3 = Vmax * (LI + L2 + 13); 

% END OF PSEUDO CODE 
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APPENDIX B 

AREA .. SUM(I,A(I)) =E= 0; 

VELOCITY(VINDX) .. VEL(VINDX) =E= VSCALE * SUM(I$(ORD(I) 

LE ORD(VINDX)), A(I)); 
POSITION .. SUM(I,VEL(I)) =E= FINALPOS * SCALEFACT; 
VLIMITP(I) .. SUM(VINDX$(ORD(VINDX) LE ORD(I)),A(I- 

(0RD(VINDX)+1))*(V0LTS(VINDX)+KBACK*VSCALE)) 
=L= VOLTLIM; 

VLIMITN(I) .. SUM(VINDX$(ORD(VINDX) LE ORD(I)), A(I- 

(0RD(VINDX)+1))*(V0LTS(VINDX)+KBACK*VSCALE)) 
=G= -VOLTLIM 

% A(I) are the current commands at time T(I) spaced equally at time DT. 

% VOLTS(VINDX) is a table of voltages representing the unit pulse response to 

% a unit output in current conunand. VOLTLIM is the voltage limit at saturation. 
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APPENDIX C 

GOALPOS .. SUM(I,A(I)*MODELAA*DT) =E=FINALPOS; 
MODEl(ILAST) .. SUM(I,-A(I)*MODELAA*MODELb/(MODELb- 

MODELa)*{EXP(-MODELa*(T(ILAST)+DT-T(I))) 

-EXP(-MODELa*(T(ILAST)-T(I))))) =E= 0.0; 
MODE2(ILAST) .. SUM(I,A(I)*MODELAA*MODELa/(MODELb- 

MODELa)*(EXP(-MODELb*(T(ILAST)+DT-T(I))) 

-EXP(-MODELb*(T(ILAST)-T(I))))) =E= 0.0; 
DERIVl(J) .. 1000.0*SUM(I,A(I)*T(I)*EXP(ZETA(J)*W(J)*T{I))* 

SIN(WD(J)*T(I))) =E= 0.0 ; 
DERIV2{J) .. 1000.0*SUM{I,A(I)*T(I)*EXP(ZETA(J)*W(J)*T(I))* 

COS(WD(J)*T(I))) =E= 0.0 ; 

% MODELAA is the mechanical gain of the system, MODELb, and MODELa 
% are the two time constants of the system in radians. One time constant is 
% associated with the L/R rise time of the motor inductance and the other is 
% the mechanical time constant of the rigid system. The A(I) are the voltages 
% which need to be determined. The T(I) are the times for each of the A(I). 
% DT is the time spacing of the outputs. W(J) are the undamped flexible 
% modes, \VD(J) are the damped flexible modes (in radians/s). 

DSI.4fi6J00.| 



- 81 - 



wo 99/45535 



PCT/US99/04910 



WHAT IS CLAIMED IS : 

1 . A method of using a position-velocity table to control a 
dynamic system, the method comprising the steps of: 

generating a position variable for the system; 

determining a velocity command for the system using the 
position-velocity table, the determining step determining the velocity command 
based on the position variable; 

shaping the velocity command in order to generate a shaped 
velocity command; and 

controlling the system based on the shaped velocity command, 

2. A method according to Claim 1, wherein the method 
controls a component of the dynamic system, the component comprising a 
head of a data storage device; and 

wherein the controlling step controls the head to move among 
various tracks of a data storage medium in the data storage device. 

3. A method according to Claim 2, wherein the generating step 
comprises comparing a preset position of the component to a measured 
position of the component ui order to determine the position variable; and 

wherein the method further comprises the step of performing 
inverse shaping on the measured position prior to comparing the measured 
position to the preset position. 

4. A method according to Claim 3, wherein the shaping step 
and the inverse shaping step reduce unwanted vibrations resulting from 
movement of the component. 
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5. A method according to Claim 3, wherein the measured 
position of the component is determined after the controlling step controls the 
component; and 

wherein the measured position of the component is fed back to 
the determining step following the controlling step. 

6. A method of generating a trajectory for inclusion in a 
position-velocity table which is used to control a dynamic system, the method 
comprising the steps of: 

generating a trajectory for the dynamic system, the trajectory 
defining system velocity in terms of system position and one or more 
additional variables; 

storing the trajectory in a position-velocity table having N 
(N>2) dimensions; and 

controlling the dynamic system in accordance with the trajectory 
stored in the position-velocity table. 

7. A method according to Claim 6, wherein the method 
controls a component of the dynamic system, the component comprising a 
head of a data storage device; and 

wherein the controlling step controls the head to move among 
various tracks of a data storage medium in the data storage device. 

8. A method according to Claim 7, wherein one of the 
variables comprises a desired movement distance of the component. 

9. A method according to Claim 7, wherein the trajectory is 
generated in real-time based on a partial fraction expansion that defines 
behavior of the dynamic system. 
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10. A method of controlling a dynamic system in accordance 
with a variation in a system variable, the method comprising the steps of: 

generating a plurality of trajectories defining system velocity in 
terms of system position, the plurality of trajectories being generated in 
accordance with at least one system variable; 

storing the plurality of trajectories in a single position-velocity 

table; 

detecting a value of the at least one system variable; and 
controlling the dynamic system in accordance with both the 
detected value of the system variable and the trajectories stored in the position- 
velocity table. 



11. A method according to Claim 10, wherein the position- 
velocity table comprises a series of trajectories corresponding to various 

15 component movement distances; and 

wherein the controlling step comprises selecting one of the 
trajectories from the position-velocity table based on the detected value of the 
system variable and controlling a component of the dynamic system in 
accordance with^the selected trajectory. 

20 

12. A method according to Claim 10. wherein the controlling 
step comprises generating a function based on the plurality of trajectories and 
the system variable, determining a single trajectory for the component based 
on the function, and controlling a component of the dynamic system based on 

25 the single trajectory. 

13. A method according to Claim 10. wherein the generating 
step comprises the steps of: 

estimating system parameters, the system parameters relating to 
30 movement of a component of the dynamic system: 
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determining whether the system parameters have varied from 
predetermined system parameters; 

modifying the trajectories based on determined system parameter 
variations: and 

storing the modified trajectories in the position-velocity table. 

14. A method of generating a trajectory for inclusion in a 
position-velocity table which is used in controlling a dynamic system, the 
method comprising the steps of: 

generating a trajectory for use in the dynamic system; 

storing the trajectory in the position-velocity table; and 

controlling the dynamic system in accordance with the trajectory 
stored in the position-velocity table; 

wherein the generating step generates the trajectory in 
accordance with a technique for reducing unwanted vibrations in the dynamic 
system. 

15. A method according to Claim 14, wherein the method 
controls a component of the dynamic system, the component comprising a 
head of a data storage device; and 

wherein the controlling step controls the head to move among 
various tracks of a data storage medium in the data storage device. 

16. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account both a system vibration limiting constraint 
and a system sensitivity constraint. 

17. A method according to Claim 16, wherein the system 
vibration limiting and sensitivity constraints reduce vibration during movement 
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of the component by less than 100%. 

18. A method according to Claim 15, wherein the technique for 
reducing unwarned vibrations of the component comprises generating the 
trajectory by taking into account one or more constraints which are a function 
of a movement distance of the component. 

19. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account a system vibration limiting constraint only. 

20. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory based on an input which has been shaped in accordance with a 
predetermined shaping function. 

21. A method according to Claim 20, wherein the input 
includes both transient portions and a steady state portion; and 

wherein only the transient portions of the input have been 
shaped in accordance with the predetermined shaping function. 

22. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory by filtering a predetermined trajectory using filters having zeros 
which are substantially near poles of the system. 

23. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account at least one of constraints relating to system 
thermal limits, system current limits, and system duty cycle. 
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24. A method according to Claim 15» wherein the technique for 
reducing unwanted movement of the component comprises the steps of: 

determining whether a trajectory excites greater than a 
predetermined level of vibrations in the system; and 

applying input shaping to the trajectory in a case that the 
trajectory excites greater than the predetermined level of vibrations. 

25. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory based on a Posicast input. 

26. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory based on a synimetric input. 

27. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating the 
trajectory based on a synmietric constraint that varies as a function of at least 
one of time and component position. 

28. A method according to Claim 15, wherein the technique for 
reducing unwanted vibrations of the component comprises generating a 
trajectory in accordance with a voltage which has been controlled by 
controlling current. 

29. A method according to any one of Claims 14 to 28, 
wherein the generating step comprises: 

identifying system parameters in real-time; and 
modifying the trajectory in real-time in accordance with the 
system parameters identified in the identifying step. 
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30. A method of controlling movement a dynamic system 
which can be expressed in terras of both rigid and flexible modes, the method 
comprising the steps of: 

generating a rigid body input for the dynamic system; 
processing the rigid body input so as to produce a processed 
input which compensates for vibrations in the flexible mode of the system; and 
applying the processed input to control the dynamic system. 

31. A method according to Claim 30, wherein the generating 
step comprises (i) creating a model of the rigid mode of the dynamic system 
based on a modal analysis, and (ii) determining the rigid body input based on 
the modal analysis. 

32. A method according to Claim 30, wherein the rigid body 
input corresponds to a fundamental limning parameter of the system, the 
fundamental limiting parameter of the system comprising a first parameter of 
the system to enter into saturation. 



33. A method according to Claim 32, wherein the processing 
step processes the rigid body input in accordance with a system vibration 
limiting constraint and a system sensitivity constraint. 

34. A method according to Claim 33, wherein the system 
vibration limiting and sensitivity constraints reduce vibration during movement 
of a component of the dynamic system by less than 100%. 

35. A method according to Claim 30, wherein the processing 
step processes the rigid body input in accordance with one or more constraints 
that are a function of a movement distance of a component of the dynamic 
system. 
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36. A meihod according to Claim 30, wherein the processing 
step processes the rigid body input in accordance with a system vibration 
limiting constraint only. 

37. A method according to Claim 30, wherein the processing 
step shapes the rigid body input using a predetermined shaping function. 

38. A method according to Claim 37, wherein the rigid body 
input includes both transient ponions and a steady state portion; and 

wherein only the transient portions of the rigid body input are 
shaped in accordance with the predetermined shaping function. 

39. A method according to Claun 30, wherein the processing 
step processes the rigid body input by filtering the input using filters having 
zeros which are substantially near poles of the system. 

40. A meihod according to Claim 30, wherein the processing 
step processes the rigid body input in accordance with at least one of 
constraints relating to system thermal limits, system current limits, and system 
duty cycle. 

41. A method according to Claim 30, wherein the processing 
step processes the rigid body input by determining a movement distance of a 
component of the dynamic system and modifying the rigid body input based on 
the movement distance. 

42. A method according to Claim 30, wherein the rigid body 
input comprises a Posicast input. 
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43. A method according to Claim 30, wherein the rigid body 
input comprises a symmetric input. 

44. A method according to Claim 30, wherein the processing 
step processes the rigid body input in accordance with a symmetric constraint 
thai varies as a function of at least one of time and position of a component of 
the dynamic system. 

45. A method according to Claim 30, wherein the rigid body 
input comprises a voltage which has been controlled by controlling current. 

46. A method according to any one of Claims 30 to 45, 
wherein the processing step comprises: 

identifying system parameters in real-time; and 
modifying the rigid body input in real-time in accordance with 
the system parameters identified in the identifying step. 

47. A method according to Claim 31, wherein the determining 
step determines the rigid body input in accordance with an insensitivity 
constraint, 

48. A method according to Claim 31, wherein the model of the 
system comprises a plurality of equations for the system; and 

wherein an insensitivity constraint for a particular system 
parameter is added to the system by taking a derivative of a system equation 
with respect to the insensitivity constraint and setting the derivative equal to 
zero. 

49. A method according to Claim 31, wherein the model of the 
system comprises a plurality of equations for the system; and 
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wherein an insensiiivity constraint for a particular system 
parameter is added to tiie system by setting a series of constraints for different 
values of the system parameter so as to limit a variation in the system 
parameter. 

50. A method according to Claim 31, wherein the rigid body 
input is determined in accordance with a feedback signal; and 

wherein the method further comprises adding a quasi-static 
correction factor to the feedback signal, the quasi-static correction factor 
correcting for a deflection in the component during movement. 

51. A method according to Claim 31, further comprising 
determining a center of mass of a component of the dynamic system; 

wherein the rigid body input is determined in accordance with a 
feedback signal based on the center of mass of the component. 

52. A method of determining plural switch times for a voltage 
input to a dynamic system having plural modes, the method comprising the 
steps of: 

generating a model of the dynamic system based on a modal 
analysis of each of the plural modes; 

determining a response of the dynamic system in terms of the 
modal analysis in the model; 

determining an expression for a contribution of each of the 
plural modes to a final location of the system based on a corresponding 
response, the contribution of each mode of the system being based on switch 
times for the voltage input; 

estimating values relating to the plural switch times; and 

calculating approximations of the values relating to the plural 
switch times based on the estimated values using the expression for the 
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contribution of each of the plural modes and the modal analysis in the model 
of the dynamic system. 

53. A method according to Claim 52, further comprising the 
step of re-calculating approximations of the values based on a previous 
approximation the values. 

54. A method according to Claim 53, wherein the re- 
calculating step is repeated a plurality of times, each lime using a re-calculated 
approximation of the values as the previous approximation of the values. 

55. A method according to Claim 52, further comprising the 
step of generating a table comprising plural switch times; 

wherein the estimatmg step comprises estimating the values 
using the table. 

56. A method according to Claim 52, further comprising the 
step of generating at least one curve corresponding to the plural switch times; 

wherein the estimating step comprises estimating the values 
using the at least one curve. 

57. A method according to Claim 52, wherein the dynamic 
system comprise a data storage device; and 

wherein the voltage inputs comprise voltage inputs to the data 
storage device. 

58. A method according to Claim 52, further comprising the 
step of performing input shaping on the voltage input after switch times 
therefor have been calculated. 
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59. A method according to Claim 52, wherein the estimating 
step is performed using a parameter estimator. 

60. A method of reducing unwanted vibrations in a dynamic 
5 system, the method comprising the steps of: 

determining whether greater dian a predetermined level of 
vibrations will be excited by a system input; and 

modifying the input to the dynamic system in a case that greater 
than the predetermined level of vibrations will be excited, where the input to 
10 the dynamic system is modified so as to reduce the level of vibrations in the 
system to less than the predetermined level of vibrations 

61. A method according to Claim 60, wherein the modifying 
step comprises using at least one of an input shaper, an inverse shaper, and a 

15 filter in order to modify the input to the dynamic system. 

62. A data storage device which uses a position-velocity table 
to control movement of a component of the data storage device, the data 
storage device comprising: 

20 a memory which stores the position-velocity table and computer- 

executable process steps; and 

a processor which executes the process steps stored in the 
memory so as (i) to generate a position variable for the component, (ii) to 
determine a velocity command for the component using the position-velocity 

25 table, the processor determining the velocity command based on the position 
variable, (iii) to shape the velocity command in order to generate a shaped 
velocity command, and (iv) to control the component to move based on the 
shaped velocity command. 
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63. A data storage device according lo Claim 62. wherein the 
component comprises a head of the data storage device; and 

wherein the processor controls the head to move among various 
tracks of a data recording medium in the data storage device. 

5 

64. A data storage device according to Claim 62, wherein, to 
generate a position variable for the component, the processor compares a 
preset position of the component to a measured position of the component; and 

wherein the processor further performs inverse shaping on the 
10 measured position prior to comparing the measured position to the preset 
position. 

65. A data storage device according to Claim 64, wherein the 
shaping and inverse shaping performed by the processor reduce unwanted 

15 vibrations resulting from movement of the component. 

66. A data storage device according to Claim 64, wherein the 
processor determines the measured position of the component after controlling 
the component; and 

20 wherein the processor uses a previously-measured position of 

the component to determine the position variable. 

67. An apparatus which generates a trajectory for inclusion in a 
position-velocity table that is used in to control a dynamic system, the 

25 apparatus comprising: 

a memory which stores computer-executable process steps and a 
position- velocity table having N (N>2) dimensions; and 

a processor which executes the process steps stored in the 
memory so as (i) to generate a trajectory for the system, the trajectory 
30 defining system velocity in terms of system position and one or more 
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additional variables, (ii) to store the trajectory in the position-velocity table, 
and (iii) to control the system in accordance with the trajectory stored in the 
position- velocity table. 

68. An apparatus according to Claim 67, wherein the apparatus 
controls a component of the dynamic system, the component comprising a 
head of a data storage device; and 

wherein the processor controls the head to move among various 
tracks of a data storage medium in the data storage device. 

69. An apparatus according to Claim 68, wherein one of the 
variables comprises a desired movement distance of the component, 

70. An apparatus according to Claim 68, wherein the processor 
generates the trajectory in real-time based on a partial fraction expansion that 
defines behavior of the dynamic system. 

71. An apparatus which controls a dynamic system in 
accordance with a variation in a system variable, the apparams comprising: 

a memory which stores a position- velocity table and computer- 
executable process steps; and 

a processor which executes the process steps stored in the 
memory so as (i) generate a plurality of trajectories defining velocity in terms 
of position, the plurality of trajectories being generated in accordance with at 
least one system variable, (ii) to store the plurality of trajectories in the 
position-velocity table, (iii) to detect a value of the at least one system 
variable, and (iv) to control the dynamic system in accordance with both the 
detected value of the system variable and the trajectories stored in the position- 
velocity table. 



- 95 - 



wo 99/45535 



PCT/US99/04910 



72. An apparatus according lo Claim 71, wherein the position- 
velocity table comprises a series of trajectories corresponding to various 
component movement distances; and 

wherein the processor controls a component of the dynamic 
system by selecting one of the trajectories from the position-velocity table 
based on the detected value of the system variable and by controlling the 
component in accordance with the selected trajectory. 

73. An apparatus according to Claim 71, wherein the processor 
controls a component of the dynamic system by generating a function based on 
the plurality of trajectories and the system variable, by determining a single 
trajectory for the component based on the function, and by controlling the 
component based on the single trajectory. 

74. An apparatus according to Claim 71, wherein the processor 
generates the plurality of trajectories by (i) estimating system parameters, the 
system parameters relating to movement of a component of the dynamic 
system, (ii) determining whether the system parameters have varied from 
predetermined system parameters, (iii) modifying the trajectories based on 
determined system parameter variations, and (iv) storing the modified 
trajectories in the position-velocity table. 

75. An apparatus for generating a trajectory for inclusion in a 
position-velocity table which is used in controlling a dynamic system, the 
apparatus comprising: 

a memory which stores the position-velocity table and computer- 
executable process steps; and 

a processor which executes the process steps stored in the 
memory so as (i) to generate a trajectory for the system, (ii) to store the 
trajectory in the position- velocity table, and (iii) to control the system in 
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accordance with the trajectory scored in the position- velocity table; 

wherein the processor generates the trajectory in accordance 
with a technique for reducing unwanted vibrations in the system. 

76. An apparanis according to Claim 75. wherein the apparaais 
controls a component of the dynamic system, the component comprising a 
head of a data storage device; and 

wherein the processor controls the head to move to among 
various tracks of a magnetic disk in the disk drive. 

77. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account both a system vibration limiting constraint 
and a system sensitivity constraint. 

78. An apparatus according to Claim 77, wherein the system 
vibration limiting and sensitivity constraints reduce vibration during movement 
of the component by less than 100%. 



20 . 79. An apparatus according to Claim 76, wherein the technique 

for reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account one or more constraints which are a function 
of a movement distance of the component. 

25 80. An apparatus according to Claim 77, wherein the technique 

for reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account a system vibration limiting constraint only. 

81. An apparatus according to Claim 76, wherein the technique 
30 for reducing unwanted vibrations of the component comprises generating the 
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trajectory based on an input which has been shaped in accordance with a 
predetermined shaping function. 

82. An apparatus according to Claim 81, wherein the input 
includes both transient portions and a steady state portion; and 

wherein only the transient ponions of the input have been 
shaped in accordance with the predetermined shaping function. 

83. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory by filtering a predetermined trajectory using filters having zeros 
which are substantially near poles of the system. 

84. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory by taking into account at least one of constraints relating to system 
thermal limits, system current limits, and system duty cycle. 

85. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted movement of the component comprises the steps of: 

determining whether a trajectory excites greater than a 
predetermined level of vibrations in the system; and 

applying input shaping to the trajectory only in a case that the 
trajectory excites greater than the predetermined level of vibrations. 

86. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory based on a Posicast input. 
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87. An apparatus according to Claim 76. wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory based on a symmetric input. 

88. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating the 
trajectory based on a symmetric constraint that varies as a function of at least 
one of time and component position. 

89. An apparatus according to Claim 76, wherein the technique 
for reducing unwanted vibrations of the component comprises generating a 
trajectory in accordance with a voltage which has been controlled by 
controlling current. 

90. An apparatus according to any one of Claims 75 to 89, 
wherein the processor generates the trajectory by (i) identifying system 
parameters in real-time, and (ii) modifying the trajectory in real-time in 
accordance with the system parameters identified by the processor. 

91. An apparatus which controls a dynamic system that can be 
expressed in terms of both rigid and flexible modes, the apparatus comprising: 

a memory which stores computer-executable process steps; and 
a processor which executes the process steps stored in the 
memory so as (i) to generate a rigid body input for the dynamic system, (ii) to 
process the rigid body input so as to produce a processed input which 
compensates for vibrations in the flexible mode of the system, and (iii) to 
apply the processed input to control the dynamic system. 

92. An apparatus according to Claim 91, wherein the processor 
generates the rigid body input by (i) creating a model of the rigid mode of the 
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dynamic system based on a modal analysis of the system, and (ii) determining 
an input to the dynamic system based on the modal analysis, 

93. An apparatus according to Claim 90, wherein the rigid 
5 body input comprises a fundamental limiting parameter of the system, the 

fundamental limiting parameter of the system corresponding to a first 
parameter in the system to enter into saturation. 

94. An apparams according to Claim 93. wherein the processor 
10 processes the rigid body input in accordance with a system vibration limiting 

constraint and a system sensitivity constraint. 

95. An apparatus according to Claim 94, wherein the system 
vibration limiting and sensitivity constraints reduce vibration during movement 

15 of the component by less than 100%. 

96. An apparatus according to Claim 91, wherein the processor 
processes the rigid body input in accordance with one or more constraints that 
are a function of a movement distance of a component of the dynamic system. 

20 

97. An apparatus according to Claim 91, wherein the processor 
processes the rigid body input in accordance with a system vibration limiting 
constraint only. 

25 98. An apparatus according to Claim 91, wherein the processor 

shapes the rigid body input using a predetermined shaping function. 

99. An apparatus according to Claim 98, wherein the rigid 
body input includes both transient portions and a steady state portion; and 
30 wherein the processor shapes only the transient portions of the 
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rigid body input in accordance with the predetermined shaping function. 

100. An apparatus according to Claim 91, wherein the 
processor processes the rigid body input by filtering the input using filters 
having zeros which are substantially near poles of the system. 

101. An apparatus according to Claim 91, wherein the 
processor processes the rigid body input in accordance with at least one of 
constraints relating to system thermal limits, system current limits, and system 
duty cycle. 

102. An apparatus according to Claim 91, wherein the 
processor processes the rigid body input by determining a movement distance 
of a component of the dynamic system and modifying the input based on the 
movement distance. 

103. An apparatus according to Claim 91, wherein the rigid 
body input comprises a Posicast input. 

104. An apparatus according to Claim 91, wherein the rigid 
body input comprises a symmetric input. 

105. An apparatus according to Claim 91, wherein the 
processor processes the rigid body input in accordance with a symmetric 
constraint that varies as a function of at least one of time and position of a 
component of the dynamic system. 

106. An apparatus according to Claim 91, wherein the 
processor processes the rigid body input based on a voltage which has been 
controlled by controlling current. 
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107. An apparatus according to any one of Claims 91 to 106, 
wherein the processor processes the rigid body input by (i) identifying system 
parameters in real-time, and (ii) modifying the input in real-time in accordance 
with the system parameters identified by the processor. 

108. An apparatus according to Claim 91, wherein the 
processor generates the rigid body input in accordance with an insensitivity 
constraint. 

109. An apparatus according to Claim 108, wherein the model 
of the system comprises a plurality of equations for the system; and 

wherein an insensitivity constraint for a particular system 
parameter is added to the system by taking a derivative of a system equation 
with respect to the insensitivity constraint and setting the derivative equal to 
zero. 

110. An apparatus according to Claim 108, wherein the model 
of the system comprises a plurality of equations for the system; and 

wherein an insensitivity constraint for a particular system 
parameter is added to the system by setting a series of constraints for different 
values of the system parameter so as to limit a variation in the system 
parameter, 

111. An apparatus according to Claim 90, wherein the 
processor generates the rigid body input in accordance with a feedback signal; 
and 

wherein the processor adds a quasi-static correction factor to the 
feedback signal, the quasi-static correction factor correcting for a deflection in 
the component during movement. 
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112. An apparatus according to Claim 90, wherein the 
processor determines a center of mass of a component of the dynamic system; 
and 

wherein the processor generates the rigid body input in 
accordance with a feedback signal based on the center of mass of the 
component. 

113. An apparams which determines plural switch times for a 
voltage input into a dynamic system having plural modes, the apparatus 
comprising: 

a memory which stores computer-executable process steps; and 
a processor which executes the process steps stored in the 
memory so as (i) to generate a model of the dynamic system in terms of a 
modal analysis each of the plural modes, (ii) to determine a response of the 
dynamic system in terms of the modal analysis in the model, (iii) to determine 
an expression for a contribution of each of the plural modes to a final location 
of the system based on a corresponding response, the contribution of each 
mode of the system being based on switch times for the voltage input, (iv) to 
estimate values corresponding to the plural switch times, and (v) to calculate 
approximations of the values corresponding to the plural switch times based on 
the estimated values using the expression for the contribution of each of the 
plural modes and the modal analysis in the model of the dynamic system. 

114. An apparatus according to Claim 113, wherein the 
processor re-calculates approximations of the values based on a previous 
approximation of the values. 

115. An apparatus according to Claim 114, wherein the 
processor re-calculates approximations of the values a plurality of times, each 
time using a re-calculated approximation of the values as the previous 
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approxinmion of the values. 

116. An apparatus according to Claim 113, wherein the 
processor generates a table comprising plural switch times; and 

5 wherein the processor estimates the values using the table. 

117. An apparatus according to Claim 113, wherein the 
processor generates at least one curve corresponding to the plural switch times; 
and 

10 wherein the processor estimates the values using the at least one 

curve. 

118. An apparatus according to Claim 113, wherein the 
dynamic system comprises a data storage device; and 

15 wherein the voltage inputs comprise voltage inputs to the data 

storage device. 

119. An apparatus according to Claim 113, further comprising 
the step of performing input shaping on the voltage input after switch times 

20 therefor have been calculated. 

120. An apparatus which reduces unwanted vibrations in a 
dynamic system, the apparatus comprising: 

a memory which stores computer-executable process steps; and 
25 a processor which executes the process steps stored in the 

memory so as (i) to determine whether greater than a predetermined level of 
vibrations will be excited by an input to the system, and (ii) to modify the 
input to the dynamic system in a case that greater than the predetermined level 
of vibrations will be excited, where the processor modifies the input to the 
30 dynamic system so as to reduce the level of vibrations in the system to less 
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than the predetermined level of vibrations. 

121. An apparatus according to Claim 120. wherein the 
processor modifies the input to the dynamic system using at least one of an 
input shaper. an inverse shaper, and a filter. 

122. A method of controlling a dynamic system in accordance 
with an input that is a function of time so as to reduce unwanted vibrations in 
the system, the method comprising the steps of: 

generating a model of the dynamic system, the model defining 
system position in terms of both tune and a system input, and the model 
constraining the system in accordance with one or more constraints relating to 
the unwanted vibrations; 

determining an input to the dynamic system which reduces the 
unwanted vibrations based on the model generated in the generating step; and 

controlling the dynamic system in accordance with the input 
determined in the determining step. 

123. A method according to Claim 122, wherein the model of 
the system comprises a partial fraction expansion of third order equations that 
define the system. 

124. A method according to Claim 123, wherein the partial 
fraction expansion equations comprise: 
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where Finalpos is the final position of a component of the dynamic system. 
Ten. corresponds to a time at which Finalpos is reached. A, a and b are based 
on the system parameters, are voltage inputs to the system, Tj are the times 
at which Vj are input, and At is a time interval at which are input. 

125. A method according to Claim 122, wherein the input 
determined in the determining step comprises the fundamental limiting 
parameter of the system, the fundamental limiting parameter corresponding to 
a first parameter in the system to enter into saturation. 

126. A method of using a current command to control a system 
having voltage as a physical limiting parameter, where the system includes a 
current controller connected to a power supply, the method comprising the 
steps of: 

inputting a current command to the system; 

shaping the current command using a unity magnitude shaper so 
that the current controller in the system goes into saturation; and 

supplying voltage to the system from the power supply via the 
current controller in saturation. 

127. A method of generating a position- velocity table for a 
dynamic system, the method comprising the steps of: 

modeling the dynamic system in terms of partial fr?ction 
expansion equations; 

integrating the partial fraction expansion equations forward in 
time so as to generate a trajectory for the dynamic system; and 

storing the trajectory for the system in the position-velocity 

table. 



- 106 - 



wo 99/45535 



PCT/US99/04910 



128. A method according to Claim 127. wherein the partial 
fraction expansion equations which model the dynamic system comprise: 

/V 

Finalpos = ^^^A/ 

i=! ' b ' a 
1=1 a - b 



where Finalpos is the final position of a component of the dynamic system. 
Tend corresponds to a time at which Finalpos is reached. A, a and b are based 
5 on the system parameters, Vj are inputs to the system, T; are the times at 
which Vj are input, and At is a time interval at which Vj are input. 



129. An apparatus according to Claim 75, wherein the position- 
velocity table comprises a non-dimensional position velocity table. 

10 

130. A method of controlling a dynamic system having one or 
more feedforward inputs, where one of the feedforward inputs corresponds to 
a fundamental limiting parameter of the system, the method comprising the 
steps of: 

15 altering a form of a feedforward input that corresponds to the 

fundamental limiting parameter of the system so as to reduce unwanted 
dynamics of the system. 

131. A method according to Claim 130, further comprising the 
20 step of determining the fundamental limiting parameter of the system by 

identifying a first parameter of the system to enter into saturation. 
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132. A method according to Claim .130, wherein the altering 
step comprises shaping the feedforward input. 

133. A method according to Claim 132, wherein the shaping is 
performed using Input Shaping^*^. 

134. A method according to Claim 132, wherein the shaping is 
performed using one or more filters. 

135. A method according to Claim 132, further comprising the 

steps of: 

identifying any nonlinear elements in the system; 
wherein the shaping is performed after any nonlinear elements 
identified in the identifying step. 

136. A method according to Claim 130, wherein the altering 
step comprises pre-saturating the feedforward input and then shaping the 
feedforward input. 



20 137. A method according to Claim 130, wherein the dynamic 

system comprises a data storage device system; and 

wherein the fundamental limiting parameter comprises voltage. 

138. A data storage device system having one or more 
25 feedforward inputs, where one of the feedforward inputs corresponds to a 

fundamental limiting parameter of the system, the system comprising: 

a memory which stores computer-executable process steps; and 
a processor which executes the process steps stored in the 
memory so as to alter a form of a feedforward input that corresponds to the. 
30 fundamental limiting parameter of the system so as to reduce unwanted 



15 
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dynamics of the system. 

139. A system according to Claim 138, wherein the processor 
executes process steps so as to determine the fundamental limiting parameter 
of the system by identifying a first parameter of the system to enter into 
saturation. 

140. A system according to Claim 138, wherein the 
feedforward input is altered by shaping the feedforward input. 

141. A system according to Claim 140, wherein the shaping is 
performed using Input Shaping^^^ 

142. A system according to Claim 140, wherein the shaping is 
performed using one or more filters. 

143. A system according to Claim 140, wherein the processor 
executes process steps so as to identify any nonlinear elements in the system; 

wherein the shaping is performed after any nonlinear elements 
identified by the processor. 

144. A system according to Claim 138, wherein the processor 
alters the feedforward input by pre-saturaiing the feedforward input and then 
shaping the feedforward input. 

145. A method of shaping an input to a dynamic system so as 
to reduce unwanted dynamics in the system, the input to the dynamic system 
comprising digital data sampled at a predetermined frequency, the method 
comprising the steps of: 
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identifying system vibrations that occur at the Nyquist frequency 
for the system, the system vibrations corresponding to a sine wave having two 
sample points per period; and 

applying a three-pulse shaper to the input, wherein first and 
5 second pulses of the three-pulse shaper are applied at the two sample points in 
a first period of the input, and a third pulse of the three-pulse shaper is applied 
at a first sample point in a second period of the input. 

146. A method of generating an input to a computer-controlled 
10 dynamic system so as to suppress vibrations therein, the dynamic system 

having a dedicated path solely for a feedforward input from a controller to 
controlled hardware, the method comprising the steps of: 

determining a frequency of vibrations to be suppressed; 
wherein, in a case that the frequency of the vibrations to be 
15 suppressed is at or below a servo rate for the dynamic system, the method 
comprises the steps of: 

executing servo calculations for the system; 
determining a servo output based on the servo 

calculations; and 

20 outputting the servo output as the input to the dynamic 

system; and 

wherein, in a case that the frequency is above the servo rate for 
the dynamic system, the method comprises the steps of: 
determining a trajectory value; 
25 shaping the trajectory; and 

outputting the shaped trajectory as the input to the 

dynamic system. 

147. A method of generating an input to a computer-controlled 
30 dynamic system so as to suppress vibrations therein, the dynamic system 
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having a path by which a feedforward input and other signals are output from 
a controller to controlled hardware, the method comprising the steps of: 

executing servo calculations for the system; 

determining a servo output based on the servo calculations; 

storing the servo output in a memory; 

determining a trajectory value for the feedforward input; 

shaping the trajectory value; and 

adding the servo output stored in the memory to the shaped 
trajectory value so as to generate the feedforward input. 

148. A graphical user interface ("GUI") which provides 
controllers for affecting operation of a data storage device, the GUI 
comprising a first controller which alters at least one of a seek time of the data 
storage device and a noise level of the data storage device. 

149. A GUI according to Claim 148, wherein the first 
controller causes progressive changes in the noise level and the seek time of 
the data storage device; and 

wherein as the first controller increases the noise level of the 
data storage device, the first controller causes the seek time of tlie data storage 
device to decrease, and as the first controller decreases the noise level of the 
data storage device, the first controller causes the seek time of the data storage 
device to increase. 

150. A GUI according to Claim 149, wherein the first 
controller comprises a sliding bar which moves along a continuum on which 
data storage device noise level and seek time vary inversely, the continuum 
including a first end comprising a high noise level/low seek time and a second 
end comprising a low noise level/high seek time. 
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151. A GUI according to Claim 150, further comprising a 
display area which displays discrete values corresponding to the noise level 
and/or the seek time of the data storage device set by the first controller. 

152. A GUI according to Claim 148. wherein the first 
controller causes the data storage device to reduce its power consumption. 

153. A GUI according to Claim 152, wherein the first 
controller causes progressive changes in the power consumption and the seek 
time of the data storage device; and 

wherein as the first controller causes the data storage device to 
decrease power consumption, the first controller causes the seek time of the 
data storage device to increase, and as the first controller causes the data 
storage device to increase power consumption, the first controller causes the 
seek time of the data storage device to decrease. 

154. A GUI according to Claim 153, wherein the first 
controller comprises a sliding bar which moves along a continuum on which 
data storage device power consumption and seek lime are inversely variable, 
the continuum including a first end comprising a high power consumption/low 
seek time and a second end comprising a low power consumption/high seek 
time. 

155. A GUI according to Claim 154. funher comprising a 
display area which displays discrete values corresponding to the seek time of 
the data storage device set by the first controller and an amount of power 
remaining in the data storage device. 

156. A GUI according to Claim 148, wherein the first 
controller comprises discrete values which are selectable to alter the noise 
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level and/or the seek time of the data storage device. 

157. A GUI according to Claim 148, further comprising a 
second controller, which is separate from the fu^t controller, the second 
controller causing the data storage device to enter into a power-saving mode, 
the power-saving mode comprising a state during which the data storage device 
reduces its power consumption by a predetermined amount. 

158. A GUI according to Claim 157, wherein the second 
controller comprises a check box. 

159. A GUI according to Claim 148, wherein settings in the 
GUI override previous settings in the data storage device. 

160. A GUI according to Claim 148, ftirther comprising a 
preview controller, the preview controller causing the data storage device to 
operate using a noise level set by the first controller. 

161. A method of controlling operation of a data storage 
device, the method comprising the steps of: 

generating a graphical user interface ("GUI"), the GUI 
providing a first controller for controlling at least one of a seek time of the 
data storage device and a noise level of the data storage device; 

operating the first controller so as to alter settings in the GUI 
for at least one of the seek time and the noise level of the data storage device; 
and 

outpuiting commands to the data storage device causing the data 
storage device to alter its operation in accordance with altered settings in the 
GUI. 
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162. A method according to Claim 161, wherein the first 
controller additionally controls power consumption by the data storage device; 
and 

wherein the operating step comprises operating the first 
controller so as to alter settings in the GUI relating to data storage device 
power consumption. 

163. A method according to Claim 161, wherein the GUI 
further comprises a second controller for controlling power consumption by 
the data storage device; and 

wherein the method further comprises, before the outputting 
step, the step of using the second controller to alter settings in the GUI 
relating to data storage device power consumption. 

164. Computer-executable process steps stored on a computer- 
readable medium, the computer-executable process steps to control operation 
of a data storage device, the computer-executable process steps comprising: 

code to generate a graphical user interface ("GUI"), the GUI 
providing a first controller for controlling at least one of a seek time of the 
data storage device and a noise level of the data storage device; 

code to operate the first controller so as to alter settings in the 
GUI for at least one of the seek time and the noise level of the data storage 
device; and 

code to output commands to the data storage device causing the 
data storage device to alter its operation in accordance with altered settings in 
the GUI. 

165. Computer-executable process steps according to Claim 
164, wherein the first controller addhionally controls power consumption by 
the data storage device; and 
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wherein the operating code operates the first controller so as to 
alter settings in the GUI relating to data storage device power consumption. 

166. Computer-executable process steps according to Claim 
165, wherein the GUI further comprises a second controller for controlling 
power consumption by the data storage device; and 

wherein the computer-executable process steps further comprise 
code to operate the second controller to alter settings in the GUI relating to 
data storage device power consumption. 

167. An apparatus for controlling operation of a data storage 
device, the apparats comprising: 

a memory which stores computer-executable process steps; and 
a processor which executes the process steps so as (i) to 
generate a graphical user interface ("GUI"), the GUI providing a first 
controller for controlling at least one of a seek time of the data storage device 
and a noise level of the data storage device, (ii) to operate the first controller 
so as to alter settings in the GUI for at least one of the seek time and the noise 
level of the data storage device, and (iii) to output commands to the data 
storage device causing the data storage device to alter its operation in 
accordance with altered settings in the GUI. 

168. An apparatus according to Claim 167, wherein the first 
controller additionally controls power consumption by the data storage device; 
and 

wherein the operating step comprises operating the first 
controller so as to alter settings in the GUI relating to data storage device 
power consumption. 
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169. An apparatus according to Claim 168, wherein the GUI 
further comprises a second controller for controlling power consumption by 
the data storage device; and 

wherein the processor executes process steps, before executing 
the outputting step, to use the second controller to alter settings in the GUI 
relating lo data storage device power consumption. 

170. A method of controlling a dynamic system using an input 
command, comprising the steps of: 

shaping the input command to saturation; 

inputting the saturated command until a first predetermined 
condition is detected; 

shaping a transition of the input command during deceleration 
from saturation until a second predetermined condition occurs; and 

following a preset trajectory until the dynamic system comes to 
widiin a predetermined proximity of its final state. 

17 L A method according to Claim 170 wherein the preset 
trajectory comprises a ciu've in a PV table. 

172, A data storage device comprising: 

a memory which stores computer-executable process steps, the 
process steps corresponding to a plurality of methods for controlling the data 
storage device with reduced vibrations; 

a manual switch for selecting process steps from the memory 
that correspond to one of the methods; and 

a processor the for executing the selected process steps so as to 
control the data storage device. 
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173. A method of generating commands for a dynamic system 
in a first parameter which maintain a limit in a second parameter, where the 
second parameter comprises a fundamental limiting parameter of the dynamic 
system, the method comprising the steps of: 
5 determining a response of the second parameter in the dynamic 

system to a unit command in the first parameter; and 

generating the command in the second parameter based on the 
response determined in the determining step. 

10 174. A method according to Claim 173, wherein the fust 

parameter is current and the second parameter is voltage; and 

wherein the dynamic system comprises a disk drive. 

175. A method according to Claim 173, wherein the response 
15 is determined by iteratively solving a set of equations for the first parameter 

knowing at least the second parameter. 

176, A method according to Claim 175, wherein the set of 
equations comprises: 

N 
1=1 

20 where A comprises amplitudes of the command in the first parameter at each 
time interval i. and N comprises a last time interval; 

/-I 

where v comprises a system velocity and Cyscaie is a constant; 

Pfinnl = E V, . 
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Pfmai comprises a final state of the system; and 



Where R comprises a pulse response of the system to the second parameter and 
^lim comprises a limit in the second parameter. 

177. A method according to Claim 176, wherein A comprises 
current, V comprises voltage, and R comprises a voltage response of the 
system. 

178. A method according to Claim 176, wherein the values of 
R(i) are determined by taking a peak value of the system response and 
sampling values of the system response at subsequent time increments. 

179. A method generating commands for a dynamic system in a 
first parameter (A) which maintain a limit in a second parameter (V), where 
the second parameter (V) comprises a ftmdamental limiting parameter of the 
dynamic system, the method comprising the steps of: 

determining a values for a command in the first parameter (A) 
at time intervals (i) based on the following relationship: 



where R comprises a pulse response of the system in the second parameter; 
and 

formulating a command over time in the first parameter (A) 
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based on the A(i) values determined in the determining step. 

180. A method according to Claim 179, wherein A comprises 
current and V comprises voltage. 

181. A method of controlling a dynamic system having 
vibrations resulting from movement, the method comprising the steps of: 

identifying transitions of an input command to the dynamic 

system; and 

shaping transitions of the input command so as to result in a 
system response to the input command with reduced vibrations. 

182. A method of controlling a system to reduce unwanted 
dynamics using commands in both first and second parameters, where the 
second parameter comprises a fundamental limiting parameter of the system, 
the method comprising: 

commanding the system in the first parameter during a first 
mode of system operation; and 

commanding the system in the second parameter during a second 
mode of system operation. 

183. A method according to Claim 182, wherein the system 
comprises a disk drive; 

wherein the first mode of operation comprises tracking 
performed by the disk drive; and 

wherein the second mode of operation comprises seeking 
performed by the disk drive. 

184. A method according to Claim 176. 178 and 179, wherein 
V,i,„ is varied in accordance with i. 
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185. A method according to Claims 173 to 179, wherein 
constraints are added for parameter slew rate limits; and 

wherein the generating step generates the command in 
accordance with the added constraints. 

186. A method of rescaling a vibration-limiting input to a 
dynamic system, the method comprising the step of: 

linearly scaling amplitudes of the vibration-limiting input to 
produce a scaled vibration-limiting input. 



- 120 - 



wo 99/45535 PCTAJS99/0491 0 

01/28 

FIG. 1 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCTAJS99/04910 






CM 




03 




«0 «0 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



05/28 



FIG. 3 




START 



S301 



GENERATE GUI 



S302 

J 
1 



INPUT SETTINGS TO GUI 



S503 



GENERATE COMMANDS BASED ON SEUINGS 



S304 



OUTPUT COMMANDS TO DISK DRIVE 




END 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



04/28 



PCT/US99/04910 



O 

o 

Q 

in 
Q 



r 



to 



o 




O 
O 

Q 



Uj Ul Ul 

5 s 2 

Uj Uj Ul 
QQQ 



-J] 



Q 
»0 



3 



Ul 



Ul 
Q 
O 



g 



7 



to 



3 

Or 



8: 



Ul 

5 



O 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



05/28 



PCT/US99/04910 




ID 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



06/28 



PCTAJS99/04910 



El 

D 



o 

z. 
o 
o 

Q 

it: 

S 




o 



o 
o 

Q 

£.5 
g 

I O 

o 



o 




Q 
O 

o 

□ 



8: 



UJ 

o 
z. 

5 



o 



CO 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



07/28 



PCTAJS99/04910 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



08/28 



PCT/US99/04910 




00 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



09/28 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCTAJS99/04910 



FIG. 1 0B 



11/28 



10 



75 




U ^74 



SENSOR 



MOTOR 



79 75 



CONTROLLER/ 
PROCESSOR 



72 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 PCT/US99/0491 0 

13/28 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



r/G. / 3 

PRIOR ART 



14/28 



SEJPOINT 



PV LOOKUP 
TABLE 



VELOCITY 



SYSTEM 

(CLOSED LOOP 
VELOCITY) 



POSITION 



FIG. 1 4 



SEJPOINT 



PV LOOKUP 
TABLE 



84 



r 



VELOCITY 



SHAPER 



V 

86 



SYSTEM 
(CLOSED LOOP 

VELOcmr) 



POSITION 



INVERSE 
SHAPER 



■87 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 PCT/US99/04910 

15/28 



FIG. 1 5 




FIG. 1 6 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 PCT/US99/0491 0 

16/28 



FIG. 1 7 



X(t) 


IS 


y(t) 







FIG 18 



x(t) 




y(t) 







SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



FIG. 1 9 



17/28 





15 









FIG20 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



20/28 






o 


■1 


1 

o 








VOLTAGE 
DLTS) 


CN 
6 






CL 











7 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 



21/28 



ii 

ii 




I 



If) 
I 



Uj 
CL 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCTAJS99/04910 



22/28 



FIG. 24 



VOLTAGE 




TIME 



FIG. 25 



"2 



o o o 



o o o 



HL 



TIME 



FIG.26 



-r — ' — I -f- o o o -f- 



TA/D 



'0 



'1 



-^2 



'N 



SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 



PCT/US99/04910 




SUBSTITUTE SHEET (RULE 26) 



PCT/US99/049I0 



25/28 



1 


1 

1 


D 








Kx 




+ 




to 




Ka 




Kv 






Ka 


(0 


Kx Kv 


i 




+ 
CO 








X 






SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



# 



wo 99/45535 



PCTAJS99/04910 



FIG. 33 



27/28 



T 1 r 



AMPLITUDE 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
TIME (sec) 



FIG. 34 



AMPLITUDE 




0.5 1 
TIME (sec) 



1.5 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45535 PCTAJS99/0491 0 

28/28 

FIG. 35 




SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



Inte lonoi Application No 

PCT/US 99/04910 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 6 G11B5/55 



According to intgrnatlonal Patent Classification (IPCI or to both national ciasaillcation and IPC 



B. FIELDS SEARCHED 



Minimum documentation searchod {ctaGSillcation system followed by ctassitication symbotei 

IPC 6 GllB 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ' 



Citation of document, with indication, whera appropriate, of thQ relevant passages 



EP 0 441 407 A (FUJITSU LTD) 
14 August 1991 (1991-08-14) 



page 5, line 16 - page 5, line 56 

-/- 



Relevant to claim No. 



1,6,10, 

14,30, 

52,60, 

62,67, 

71,75. 

91,113, 

120.122, 

126,127, 

130.138, 

145-148, 

161,164, 

167,170. 

172,173. 

179,182. 

186 



)( Further documents are listed in the continuation of box C. 



ID 



Patent family members are listed in annex. 



" Special categories of cited documents : 

"A" document detining the general state of the art which is not 
considered to be of particular relevance 

"E" earlier document but published on or after the international 
liiing date 

"L" document which may throw doubts on priority ctaim(s) or 
which is cited to establish the publicalion date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use. exhibition or 
other means 

"P" document published prior to the intsmational filing date but 
taier than the priority date claimed 



T' later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed inventiw 

cannot be considered lo involve an Inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious lo a person skilled 
in the art. 

document mennber of the same patent family 



Date of the actual completion ol the international search 



22 July 1999 



Date of mailing of the international search report 



29/07/1999 



Name and mailing address of the ISA 

European Patent Offico, P.B. 58 18 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel, (+31-70) 340-2040, Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Authorized officer 



Ressenaar, J-P 



Foim PCT/I&A/210 (second shoei) (July i992) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



Inte onal Application No 

PCT/US 99/04910 



C.(Continuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 




Citation ol documsnt, wilh in{liC3tiori,wh6fo spproprists, of tho rolBVdni pdss^QOS 


R6l8vant to claim No. 


X 


us 5 696 647 A (PHAN DUC T ET AL) 


1.6,10, 






1 A OA 

14 , 30 , 






52.60, 






62.67, 






71.75, 






91,113, 






120,122, 






126,127, 






130,138, 






145-147 - 




column 6, line 59 - column 8, line 57 
— 




A 


EP 0 543 654 A (FUJITSU LTD) 


1.6.10, 




26 May 1993 (1993-05-26) 


14,30, 






b<;,oO, 






62,67, 






71,75, 






91,113, 






120,122, 






126,127, 






130,138, 






145-147 




page 4, line 2 - page 5, line 41 






page 7, line 33 - page 13, line 25 

— 




A 


EP 0 308 062 A (IBM) 


1,6,10, 




ZZ March 1989 (1989-03-22) 


14,30, 






52,60, 






62,67, 






71.75, 






91.113. 






120,122, 






126,127, 






130,138, 






145-147 




column 5, line 33 - column 7, line 38 




A 


— 

GB 2 303 732 A (SAMSUNG ELECTRONICS CO 


1,6,10, 




LTD) 26 February 1997 (1997-02-26) 


14.30, 






52.50, 






62,67, 






71.75, 






91,113, 






120,122, 






126,127, 






130,138, 






145-147 




page 13, line 38 - page 15, line 21 





Fbtm PCT/ISA/21 0 (continuation of second sheet) (July 1992} 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

Infonnatien on patent family members 



Inlft. .lonal Application No 

PCT/US 99/04910 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member{s) 



Publication 
date 



EP 


0441407 


A 


14- 


-08- 


■1991 


JP 


2657561 B 


24-09-1997 














JP 

u r 


32'i3509 A 


17-in-lQQl 

1/ iV l77i 














CA 


2036024 A 


l\J VO 17 7 1 














DE 


69115944 0 
















DE 


69115944 T 


17 t/7 177U 














KR 


9309998 B 


IJ iU 1770 














US 


5151639 A 


C7 U7 177£ 


us 


5696647 


A 


09- 


-12- 


•1997 


US 


5726825 A 


10-03-1998 














US 


5760992 A 


02-06-1998 














US 




1 o.nc^-IQQA^ 

KiO l77Q 


EP 


0543654 


A 


26- 


-05- 


•1993 


JP 


^ / *J\J / LiJ O 


Utl UH l770 














JP 




11 UO 1!770 














JP 




Ud UH i770 














JP 


5158543 A 


25-06-1993 














DE 


69227434 D 


03-12-1998 














DE 


692274'54 T 
















KR 


9613213 B 


02-10-1996 














US 


5469414 A 


21-11-1995 


EP 


0308062 


A 


22- 


■03- 


•1989 


US 


4816941 A 


28-03-1989 














CA 


1317372 A 


04-05-1993 














DE 


3875439 A 


26-11-1992 














JP 


1169785 A 


05-07-1989 














JP 


2111552 C 


21-11-1996 














JP 


8028072 B 


21-03-1996 


GB 


2303732 


A 


26- 


■02- 


•1997 


DE 


19615964 A 


30-01-1997 



R>rm PCT/rs/Vgio (patent family annex) (July 1992) 



I 



